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A COMPUTER PROGRAM FOR COMPOSING COMPRESSOR BLADING FROM 
SIMULATED CIRCULAR-ARC ELEMENTS ON CONICAL SURFACES 

by James E. Crouse, David C. Janetzke, 
and Richard E. Schwirian 

Lewis Research Center 
SUMMARY 

In axial flow compressors, the design blade elements lie on conical surfaces which 
approximate the actual stream flow surfaces. A blade-element-layout method is de- 
veloped which preserves the constant-angle change characteristic of the circular-arc 
profile. More specifically, the mean camber line and the suction and pressure surface 
lines of a blade element are lines with a constant rate of angle change with path distance 
on a specified conical surface. The layout method developed in this report and incor- 
porated in a computer program has the capability of handling a multiple-circular-arc 
blade element. A .complete blade is composed by stacking design blade elements on a 
line which may be tilted in the tangential and axial directions to minimize blade stresses. 

Blade surface coordinates for plane sections through the blade are computed for use 
in the mechanical design and fabrication of the blade. The area, center of area, and 
moments of inertia for each blade section are computed for use in stress analysis of the 
blade. 


INTRODUCTION 

In present-day aircraft compressors, the annulus area converges sufficiently through 
a blade row that the meridional (radial-axial plane) streamlines near the hub and/or tip 
have significant slopes. Compressor blading is usually constructed from blade elements 
which are designed to turn the flow on the meridional streamlines. Each individual 
blade element is generally assumed to lie on a conical surface representation of the axi- 
symmetric streamsurface through a blade row (see fig. 1). However, the layout of a 
blade element on a cone cannot retain all the properties of a conventional blade shape 
(e. g. , double circular arc). The blade -element-layout problem is to preserve the de- 




sirable properties of a conventional blade shape. There is no standard method of simu- 
lating conventional blade shapes on a cone. Commonly used methods are (1) layout of a 
reference blade element directly on the conical surface and (2) projection of a reference 
blade element laid out on a plane or a cylindrical surface to the conical surface. With 
low streamline slopes, each of these methods gives essentially the same shape on a 
cone. However, with large streamline slopes, these methods can give significantly dif- 
ferent blade shapes on a cone. 

Blade surface velocities and pressures are influenced by several interacting forces, 
but probably the most direct factor controlling local flow on a blade element is the 
rate of surface-angle change (surface curvature) along the streamline. Then, perhaps, 
the most fundamental method of simulating a desirable blade element is to retain the rate 
of surface-angle change. This approach is taken in this report to develop a computerized 
method for simulating circular-arc-type blade elements. 

The design of high-speed compressors has made wide use of blade rows composed of 
double -circular -arc (DC A) blade elements. A DCA blade element consists of one circu- 
lar arc forming the suction surface and another forming the pressure surface. This type 
of blade element has performed very successfully, and extensive data from both two- 
and three-dimensional cascades has been incorporated into blade design procedures 
(ref. 1). 

More recently, the need to control shock loss and throat area in the blade passages 
of transonic compressors has led to the use of multiple-circular-arc (MCA) blade ele- 
ments. An MCA blade element consists of two circular arcs forming the suction surface 
and two others forming the pressure surface. This type of blade permits additional con- 
trol of the chordwise turning (loading) distribution and aids in controlling the shock loss 
in blade passages with supersonic flow (refs. 2 to 5). 
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The main part of this report is a detailed development of a layout method which 
simulates an MCA blade element on a conical surface. The developed blade-element- 
layout method preserves the constant rates of angle change of the centerline and the sur- 
faces. Following the layout-method development, a step-by-step procedure for compos- 
ing a complete blade by stacking the blade elements is given. The layout method and the 
stacking procedure are incorporated in a computer program to calculate the coordinates, 
areas, and other related properties of the blade cross-sections. This computer program 
eliminates the lengthy graphical procedures previously used in the mechanical design of 
a compressor blade. 


COMPARISON OF SOME LAYOUT METHODS 

In order to illustrate the potential effect of layout method on the rate of angle change 
of a blade-element centerline, a comparison of five layout methods is presented. The 
differences are best shown with an example of a DCA blade element at the hub of a com- 
pressor. The blade parameters selected and held constant on the hub cone are the fol- 


lowing: 

Streamline slope in the meridional (r-z) plane, cx, deg 45 

Ratio of blade-section outlet radius (trailing edge) to inlet 

radius (leading edge) , r /r^ 1.4 

Leading-edge blade angle, k^, deg 45 

Trailing-edge blade angle, k q , deg 0 


These parameters and other nomenclature for the layout on a cone are shown in figure 2. 
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Blade elements which have circular -arc centerlines on a plane were laid on the cone 
by using the following layout methods: (1) a constant rate of change of local blade angle 
on the cone with distance (constant d/c/ds), (2) a circular-arc element laid on a cone, 

(3) a circular-arc element laid on a plane perpendicular to the stacking axis and pro- 
jected to the cone by lines parallel to the radial stacking axis, (4) a circular-arc element 
laid on the cylinder of blade-element outlet radius and projected to the cone by lines par- 
allel to the radial stacking axis, and (5) a circular-arc element laid on the cylinder of 
blade-element outlet radius and projected to the cone by radial lines from the axis of ro- 
tation. 

The rates of change of local centerline blade angle with distance along the blade - 
element centerline on the cone d/c/ds were computed for each of the layout methods. 

The results are compared in figure 3. Each curve has the same /c change from inlet to 
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Axial distance along blade from inlet to outlet, z/z 0 

Figure 3. - Comparison of circular-arc-layout methods. 


outlet; but the s distance varied slightly to match the specified radius change and cone 
angle conditions. The values of d/c/ds are negative since k decreases with s. 

The line of constant d/c/ds in figure 3 is from the first layout method. With all of 
the other layout methods, k changes the most rapidly at the leading edge of the blade. 
The percentage increases of d/c/ds for the other layout methods at the blade leading 
edge, as compared with the constant d/c/ds method (method 1) are 23 percent for the 
geometric layout (method 2), 35 percent for the parallel projection methods (methods 
3 and 4), and 127 percent for the radial projection method (method 5). Figure 3 shows 
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that the blade-element-layout method can have a significant effect on the d/c/ds prop- 
erties of a blade airfoil. If an important blade shape property, such as d/c/ds, is 
changed significantly by the layout method, similar changes in blade-element perform- 
ance may also be expected. 


DEVELOPMENT OF EQUATIONS FOR BLADE-ELEMENT LAYOUT 

The layout of blade elements is one of the latter steps of a compressor design. In 
the steps preceding the layout, the selections of values for the blade-element properties 
are made . For the purpose of this report, the following values are presumed to have 
been established: 

(1) Radial distance from the axis of rotation to the leading-edge center, r^ c 

(2) Radial distance from the axis of rotation to the trailing-edge center, r QC 

(3) Thickness at the leading edge, t^ 

(4) Thickness at the maximum thickness point, t 

(5) Thickness at the trailing edge, t 

(6) Angle of the centerline at the leading edge, /c^ c 

(7) Angle of the centerline at the transition point, k ^ 

(8) Angle of the centerline at the trailing edge, /c QC 

(9) Axial distance from the leading-edge center to the maximum thickness point on 

the centerline, z mc - z^ c 

(10) Axial distance from the leading-edge center to the transition point on the 
centerline, z^ c - z- c 

(11) Axial distance from the leading-edge center to the trailing-edge center, 

z - z. 
oc 1C 

These parameters and some of the nomenclature used to describe the blade elements are 
shown in figures 4 and 5. 

In the following development the constant d/c/ds property of the MCA blade element 
is preserved in the layout onto a conical surface. The centerline, pressure surface, 
and suction surface are each formed by two segments, an inlet segment and an outlet 
segment. Each segment has its own constant d/c/ds value which, generally, is differ- 
ent from that of any other segment. The development and the forms of the equations 
used in a computerized MCA blade -element-layout method are given in the following 
sections . 
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Coordinate System for Blade Element 

The most convenient coordinate system for describing a blade element on a cone is 
the R-e system shown in figure 2. Since a cone is a single curved surface which can 
be unwrapped on a plane, the following development can be considered to be carried out 
on a plane with the R-e coordinate system of a cone. In the R-e system, R is the 
length of a ray from the cone vertex to an arbitrary point, and e is the angle from a 
reference ray to a ray passing through the arbitrary point. 
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Mathematical Description of Constant-Turning-Rate Segment 


The blade angle k is the angle between the ray R and a tangent to the blade 
centerline or surface path s. For a fixed turning rate, k decreases at a constant rate 
C as s is increased; that is, 


d k _ 
ds 


or 


ds = - 


d/c 

C 


(i) 


From figure 4 note that 


dR = cos k ds 


Rde = sin k ds 


Substituting ds from equation (1) into equations (2) and (3) gives 


dR = 


cos k d/c 


de = - sin - d/c 
RC 


( 2 ) 

(3) 

(4) 

(5) 


Equation (4) integrated is 


R - R, = — (sin Ki - sin k) (6) 

1 C 1 

where the subscript 1 refers to a point where R and k are known. (All symbols are 
defined in appendix A.) By regrouping the terms in equation (6), a particular constant £ 
is formed. 

£ = RC + sin k = RjC + sin (7) 
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Solving equation (7) for R gives 


R _ £ - sin k 
C 

Equation (8) gives R as a function of k on a segment with known constants, C and 
The differential equation for e is obtained by the substitution of equation (8) into 
equation (5). 


de = — s * n - d K (9) 

sin k - C 

However, if C = 0, k is a constant, and the following differential equation for e ap- 
plies: 


de = tan k — 
R 


( 10 ) 


In general, e is given by 


e= € 1 + i(K,K 1 ,^,R,R 1 ) (11) 

where f(x, x^, ?, R, R^) is the integral of equation (9) if C 4- 0, or equation (10) if C = 0. 
The integral of equation (9) has three solutions dependent on the value of £. Details of 
the solutions for f(x,x^, £, R, Rj) are presented in appendix B. 


Definition of Blade-Element Centerline 

In this blade-element-layout procedure, it is first necessary to establish the blade - 
element centerline. Desired blade properties (e. g. , k j, k k q , C., C ) are generally re- 
lated to the centerline, and the blade thickness is applied to the centerline. 

In this development, the blade-element centerline is composed of two constant 
dx/ds segments, an inlet segment and an outlet segment. These segments are tangent 
at a point called the transition point (see fig. 4). 

In order to determine the R-e coordinates of the centerline segments, it is first 
necessary to calculate the cone half-angle a. From the input data. 
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( 12 ) 


a 



The R coordinate of the leading-edge center is given by 


R; = 


r ic 


1C sin a 


(13) 


Note that oi = 0 cannot be used in equation (13). However, a separate set of equations 
for this special case is not warranted. A sufficiently equivalent blade element can be 
calculated by using a small cone half-angle (a = 0. 1°). 

The R coordinates of other points specified on the centerline are determined by 
equation (14) 


R = 


R ic + 


z - z ic 

COS Oi 


(14) 


After the R coordinates of the leading- edge point, the transition point, and the 
trailing-edge point on the centerline are determined, the turning constants for both 
centerline segments can be calculated. Since the blade angles at the endpoints of these 
segments are given in the input, the turning-rate constant C for a segment is obtained 
from a rearrangement of equation (6) 


sin k 1 - sin k 

C = (15) 

R - Rj 

It should be noted that the value of C for the centerline is calculated rather than speci- 
fied. The reason for this is that small errors in C and in k for small values of C 
will produce large errors in R^ c and R oc * Thus, the relative axial locations of the 
segment endpoints are specified instead. 

For convenience, the angular coordinates of a blade element are referenced from 
the leading-edge center (i. e. , e. = 0). The angular coordinates of the other endpoints 
of the two segments are calculated by equation (11). 
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Determination of Blade-Element Surfaces 


The blade-element surface curves are also composed of two constant cte/ds seg- 
ments. The pressure surface and the suction surface each have an inlet segment and an 
outlet segment which are joined at a transition point. These surface curves must satisfy 
the tangency requirement at the transition point and the thickness specifications. The 
thickness is specified at three points: the leading edge, the trailing edge, and the maxi- 
mum thickness point. 

The constants for each surface segment are determined from two points on the seg- 
ment and the slope at one of these points. The general equations needed and the methods 
used for calculating the coordinates of these points, the surface slopes at these points, 
and the resulting constants for each segment are given below. Specific applications of 
these equations are given in appendix B. 

The initial points for establishing the surface curves are calculated by applying the 
thickness specifications at three points on the centerline: the leading edge, the maxi- 
mum thickness point, and the trailing edge. On a plane surface, thickness is generally 
measured along a line perpendicular to the blade centerline. On the conical surface, the 
thickness path is described by a constant angle K n path which is normal to the center - 
line, where 


K n = K c ± g (16) 

(see fig. 5). The plus sign in equation (16) gives the path direction to the suction sur- 
face, and the minus sign gives the path direction to the pressure surface. 

The differential equations for this slightly curved thickness path are 

clR = cos K n d-> (17) 


and 


Rde = sin «: n d (18) 

where A is the distance from the centerline, as shown in figure 5. Equation (17) inte- 
grates to 


R - R = A cos K 
c n 


(19) 
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Elimination of d<J by combining equations (17) and (18) gives 


de = tan K n 


dR 

R 


( 20 ) 


Equation (20) integrates to 


e - £ c = tan K n (21) 

For the special case of k q = 0, equation (21) becomes indeterminate because K n = ± tt/2, 
and equation (19) yields R = R c> Since R is constant for this case, equation (18) can 
be directly integrated to give 


€ - 



( 22 ) 


where the plus sign is for K n = tt/2 and the minus sign is for k = -n/2. 

The coordinates of the leading edge, the maximum thickness point, and the trailing 
edge on the suction surface and the pressure surface are determined by equations (19) and 
either (21) or (22). The necessary R c values for these equations are determined from 
the input by equations (13) and (14). The k c values for the leading and trailing edges 
are given in the input. The k q value for the maximum thickness point is determined 
from the rearrangement of equation (7), which gives 

k = sin _1 (? - RC) (23) 


where £ and C are constants of the segment containing the point. 

The k angle at the maximum thickness point on either the suction or pressure sur- 
face is equal to the k angle at the maximum thickness point on the centerline, or 


K mc K ms K mp K m 


(24) 


This angle k , along with the coordinates of the maximum thickness point and either 
the leading-edge point or the trailing-edge point, provides sufficient conditions for es- 
tablishing the surface curve for the segment containing the maximum thickness point. 
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Transition line 



(a) Case 1: coincident maximum thickness and transition points. 



(b) Case 2: maximum thickness behind transition point. 



(c) Case 3: Maximum thickness ahead of transition point 


Figure 6. - Locations of maximum thickness with respect to 
transition point. 
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To permit design flexibility, the maximum thickness point can be located on either 
segment, inlet or outlet, or at the transition point. These three cases, as shown in fig- 
ure 6, are 

(1) Maximum thickness at the transition point 

(2) Maximum thickness on the outlet segment behind the transition point 

(3) Maximum thickness on the inlet segment ahead of the transition point 

In establishing the surface equations, the calculations begin on the segment contain- 
ing the maximum thickness. On this segment, two points and a slope are known for both 
the pressure surface and the suction surface. Use of these known surface conditions in 
equations (7), (11), and (15) gives three equations with three unknowns: £, C, and k . 
Elimination of £ and C leaves an equation with one unknown, k . However, the com- 
plexity of this equation makes it difficult to solve explicitly. So an iterative method is 
used to solve for £, C, and k. This iterative method consists of estimating k and 
checking the resulting e-coordinate with the known e-coordinate. 

The next step in establishing the surface equations is the calculation of the transition 
point on both the pressure surface and the suction surface. This calculation involves 
finding the intersection of the surface curves with the thickness path of the transition 
point. Use of the known conditions in equations (11), (21) or (22), and (23) gives three 
equations with three unknowns: k , R, and e. Again, the complexity of equation (11) 
makes it difficult to solve for the unknowns explicitly. An iterative method is used to 
solve for the unknowns. This iterative method consists of estimating R, then comparing 
a calculated R^ c with the known R tc The calculated R^. c is determined by a re- 
arrangement of equation (21) 


R tc (R) = 


R 


exp [tan K tc (€ tc - €)] 


(25) 


This step does not apply to case 1, where the maximum thickness and the transition 
points coincide. 

The final step in establishing the surface equations is to obtain the unknowns C, 
and k for the surfaces of the remaining segment. Two points, the transition point and 
either the trailing-edge or the leading-edge point, and the angle at the transition point 
on the pressure surface and the suction surface of the remaining segment are now 
known. The iterative method used in the first step is used in this step to calculate the 
final unknowns. 
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DESCRIPTION OF COMPLETE BLADE 


The complete blade is described from a selected number of blade cross-sections. 
These cross-sections, hereinafter called blade sections, lie on planes perpendicular to 
a radial line. The blade-section surface coordinates are obtained by stacking the blade 
elements (which lie on conical streams urfaces) in a suitable manner and fairing between 
them. A primary objective in the stacking process is to minimize blade stresses. This 
is accomplished by allowing the straight stacking line to be leaned (from a true radial 
line) at prescribed angles in both the tangential and axial directions. 

The blade-element stacking procedure requires an iterative positioning of the blade 
elements until the centers of area of the blade sections are coincident with the stacking 
line within a given tolerance. The specific steps used in the stacking procedure are the 
following: 

(1) Initial positioning of blade elements along the stacking line. The intersections of 
the stacking line with the conical streamsurfaces of the blade elements are called the 
blade-element stacking points. For the first iteration the stacking points are located at 
the centers of area of the blade elements. 

(2) Calculation of stacking points relative to common reference. The coordinates of 
the blade-element stacking points are translated into cylindrical coordinates and refer- 
enced from the center of the leading edge of the hub blade element. 

(3) Calculation of blade-element coordinates. Blade-element surface coordinates in 
a Cartesian coordinate system (x-y-z) are calculated at specific z values for fairing 
convenience. 

(4) Calculation of blade-section coordinates. Blade sections lying on planes through 
each blade-element stacking point are defined. The surface coordinates of a blade sec- 
tion are obtained from the intersections of the plane of the blade section and the 

z -fairings of the blade -element surface coordinates. 

(5) Calculation of centers of area of blade sections. The center of area for each 
blade section is found by integrating over the area defined by the blade-section coordi- 
nates . 

(6) Calculation of new blade-element stacking points. A new stacking point for each 
blade element is obtained from the intersection of a line faired through the centers of 
area of the blade sections and the conical streamsurface of each blade element. If the 
new stacking points are sufficiently close to the old stacking points, the stacking proce- 
dure is considered to be converged or finished. If not, the procedure is repeated start- 
ing at step 2 using the new stacking points. 

The blade-element stacking procedure, including these steps, is described in detail 
following a description of the coordinate systems used. 
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Coordinate Systems for Complete Blade 


In addition to the R-e coordinate system for the blade elements, two other coordi- 
nate systems are used in the stacking procedure for describing the complete blade. A 
cylindrical coordinate system (r-0-z) is used for describing the stacking line and alining 
the stacking points of the blade elements along the stacking line (see fig. 7). A Cartesian 
coordinate system (x-y-z) is used for obtaining plane sections of the complete blade (see 
fig. 8). The z-axis is common to both systems and lies along the machine axis of rota- 
tion. The direction of the z -axis is defined as positive from blade inlet toward blade out- 
let. The origin, z = 0, is defined by the axial location of the center of the hub-blade - 
element, leading-edge radius. 

The orientation of the cylindrical coordinate system is shown in figure 7. The an- 
gular coordinate 9 is measured from the r-z plane which contains the hub -blade - 
element, leading -edge center. The positive 9 -direction is from the blade pressure 
(lower) surface toward the blade suction (upper) surface. 





Figure 8. - Cartesian coordinate system for blade. 
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The orientation of the Cartesian coordinate system is shown in figure 8. The x-axis 
is parallel to the radial line which passes through the hub-blade -element stacking point. 
The positive x-direction is from hub to tip. The positive y-direction is from the blade 
pressure surface toward the blade suction surface. 


Stacking Procedure 


The objective of the stacking procedure is to position each blade element such that 
the centers of area of all blade sections lie on the stacking line. The steps in the iter- 
ative procedure are as follows: 

Initial positioning of blade elements along stacking line . - The first step in the 
stacking procedure is the initial positioning of the blade element along the stacking line. 
A good first approximation to the desired stacking of the blade elements is obtained by 
alining the centers of area of the blade elements along the stacking line. A sufficiently 
accurate calculation for these centers is given by the following equations: 


R 


sp 


e 


sp 


/ 


R dA 


f dA 

/ £dA 

f dA 


( 26 ) 


( 27 ) 


where 
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( 30 ) 



The integrals in equations (26) to (30) are evaluated by numerical integration techniques 

since the functions e„(R) and e (R) are very difficult to integrate. 

s p 

Calculation of stacking points relative to a common reference . - The second step in 
the stacking procedure is the calculation of the cylindrical coordinates of the blade- 
element stacking points relative to a common reference. For convenience, the reference 
for the 9-z coordinates is the center of the leading -edge radius of the hub blade ele- 
ment. 

The simplest and perhaps most commonly used stacking line is a radial line. How- 
ever, when blade stresses are high, the maximum blade stress can be lowered by leaning 
the stacking line slightly to introduce a centrifugal force bending moment to counterbal- 
ance the aerodynamic blading moment. In this report the stacking line is a straight line 
which can be leaned in both the 9 -direction and the z -direction from a radial line at the 
hub -blade -element stacking point. The lean angle 77 is positive in the positive 
6-direction, and the lean angle \ is positive in the positive z -direction (see fig. 7). 

From geometric considerations, it can be shown that the blade -element stacking 
point location on the stacking line is given by 


sp 


= z 


sp, h + 


(r 


sp 


r sp, h> tan x 


(31) 


0 sp = B 


sp, h + 


(32) 


where 


r sp = R sp sln “ 

(33) 

0 - £ sp. h 

S P> k sin Qfj^ 

(34) 



(35) 
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The h subscript refers to the hub-blade-element values. The lean angles, 77 and A, 
are input information for the computer program and, therefore, are presumed to have 
been calculated or estimated. 

Calculation of blade-element coordinates . - The third step in the stacking procedure 
is calculation of the x-y-z coordinates of the blade elements. The general conversion 
equations for calculating x-y-z coordinates from the R-e coordinates are 


x = R sin 


a cos/— - — + 9. 

\sin a 1 

y = R sin a sin/— - — + . ) 

Una 1C sp ’ h / 


J sp '* v sp 


%,h) 

(36) 

0 sp,h) 

(37) 

a 

(38) 


where the z value of the stacking point is z oi _ and the cylindrical coordinate of the 

sp 

center of the leading -edge radius is 


^ic ^sp ” ^sp 


e. ) = 9 , + 6 - 

ic' sp, h 


sp 
sin a 


(39) 


Since the R-e corrdinates of the leading -edge, transition, maximum thickness, and 
trailing -edge points on the blade -element surfaces have been calculated previously, the 
x-y-z coordinates of these particular points can be calculated directly with equations 
(36) to (38). 

The blade surface curve fits are most conveniently made at constant z values. 
However, before particular values of z are determined, the maximum z range for the 
complete blade is found. The minimum z is found by searching the leading-edge coor- 
dinates of both surfaces of all blade elements. The maximum z is found by the same 
type of search on the trailing-edge surface coordinates. Then, equally spaced z values 
are calculated to cover the complete z range of the blade. 

Before the surface x and y coordinates can be found, it is necessary to calculate 
the surface R and e values at the prescribed z values. For a given element the R 
coordinate is given by equation (14). Before the e-coordinate can be found, the surface 
tangent angle k is calculated by equation (23). The e -coordinate then is given by equa- 
tion ( 11 ) when the known values at the transition point are used for reference values. 
Finally, the x-y coordinates are calculated by the general equations (36) and (37), for 
each z value on each blade element to complete the information needed for the curve 
fits across the blade elements. 
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Calc ulation of blade-section coordinates . - The fourth step in the stacking procedure 
is interpolation of the blade-element surface coordinates to define blade sections. Each 
blade section has a constant x value, so the y-z surface coordinates define the blade- 
section profile. The x values used in the program are the x-coordinates of the pre- 
viously calculated stacking points of the blade elements. 

A second-order Lagrangian interpolation technique is used to calculate y^ for a 
given x on each of the surface lines of equal z values. The blade-element coordinates, 

?l,(p,s)’ y l,<p,s)’ z ]]- l?2, (p, s)’ y 2, (p,s)> z j]' and [ x 3,(p,s), y 3,(p, s )> z j} are c ° nse “>- 
tive points along a z -value line. The x-dimension falls within this group of points. The 

equation for y^ g ^ is 


y (p, s) = y l, <p, s) W l + y 2, (p, s) w 2 + y 3, (p. s) W 3 


(40) 


where 


and 


W i = 



W 3 = 



x - X 2, <p, sj 1 

_ x - x 3,(p,s)_ 


[ X l, (p,s) " x 2, (p,s) 

|[ X UP,s)- X 

3, (P, -)] 



x • x i, (p, J 

X_X 3,(p,s). 


[ X 2, (p, s) ” x l, (p, s)] 

x 2,(p,s) " x 

3, (P, s)J 



" x - x i, <p, J 

x ' x 2, (p,s). 


[ X 3, (p, s) ' x l, (p, s) 

x 3,(p,s) " x 

2, (P, s)] 


(41) 


(42) 


(43) 


The coordinates y and y are calculated for each z value. 

P ® 

Calculation of centers of area of blade sections. - The fifth step in the stacking 
procedure is the calculation of the center of area of each blade section. The coordinates 
of the center of area are determined by dividing the area moments of the blade section 
by the area of the blade section. Both the area moments and the area of the blade sec- 
tion are determined by numerical integration. 

The equations for the area and the area moments of a blade section are as follows: 



(44) 
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(45) 


y 




z 



z dy dz 



(46) 


Calculation of new blade -element stacking points . - The sixth step in the stacking 
procedure is the calculation of the new blade -element stacking points. A new stacking 
point for each blade element is calculated by curve -fitting the center -of -area coordinates 
of the blade sections and finding the intersections of the curve -fit with the conic stream - 
surface of each element. The first approximation for the new y gp and z gp coordinates 
of a new stacking point is made by interpolating the center -of -area coordinates at the old 
x_ . Then, using the y_„ and z CKi approximations, an approximate x__ is calculated 

op op op op 

from the following equations: 


r _ r + Z sp ~ Z sp, old 

sp- sp, old tana 


(47) 


x 


sp 



(48) 


The approximate x g p is then used to interpolate the center -of -area coordinates for the 

new y_ and z_„. A new x__ is calculated by using the new y c „ and z in equa- 
sp sp sp sp sp 

tions (47) and (48). 

To determine whether repositioning of the blade elements is necessary, the absolute 

differences between the old and new y and z coordinates are summed in the 

sp sp 

manner of the following equation, and the sum is compared to the specified tolerance 
limit given in the input. The equation for summing the differences is 


® ^ (| y sp, new y sp, oldj + | z sp, new z sp, old|) ^ 

i=l 

where n is the number of blade elements. If S is within the specified tolerance limit, 
the stacking procedure is considered to be converged or finished. 
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If repositioning is necessary, the new stacking point coordinates, x , y , and 

sp sp 

z g p, are used to calculate the cylindrical coordinates of the new blade -element stacking 
points: 



0ic )old 


- 6 + tan 



-1 y sp 


(50) 


R = R. + 

sp, new ic 


z - z. 
sp ic 

cos a 


(51) 


r 


sp 



(52) 


The cylindrical coordinates of the new stacking points are used in the second step of the 
stacking procedure to begin another iteration. 


FINAL CALCULATIONS AND OUTPUTS 

The final calculations and outputs of the computer program are primarily intended 
for use in the mechanical design and fabrication of a compressor blade. However, the 
calculated parameters and coordinates of the blade elements may be of interest in an 
analysis of the aerodynamic design. For this purpose, the parameters and coordinates 
of the blade elements are printed out. 

The blade-element parameters printed out are 

(1) Cone half -angle , a 

(2) Blade angle at the maximum thickness, k 

(3) Centerline blade angles at the leading edge, the transition point, and the trailing 

edge, K ic , K tc , and k qc 

(4) Pressure surface blade angles at the leading edge, the transition point, and the 

trailing edge, « ip , « tp , and K Qp 

(5) Suction surface blade angles at the leading edge, the transition point, and the 

trailing edge, K is , K ts , and n QS 

(6) Inlet and outlet segment turning rates for the centerline, C ic and C QC 

(7) Inlet and outlet segment turning rates for the pressure surface, C^ p and C Qp 

(8) Inlet and outlet segment turning rates for the suction surface, C. and C 

is qs 

The blade -element coordinates printed out define the surface profile and locate par- 
ticular points of the blade elements. The coordinates which define the surface profiles 
of the blade elements are given as x and y for the suction surface and the pressure 
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surface at a z value for each element. These are the blade-element coordinates which 
are curve -fit to obtain the blade -section coordinates. The x-y-z coordinates of parti- 
cular points are given at the leading-edge point, maximum thickness point, transition 
point, and trailing -edge point on the suction surface, pressure surface, and centerline 
for each element. 

The blade -section coordinates are the primary output of the computer program. 

The locations, or x values, of the blade sections are determined in the stacking proce- 
dure or can be specified in the input. The blade sections are described in two separate 
sets of coordinates. One set is called the unrotated coordinates and uses the x-y-z 
coordinate system of the stacking procedure (see fig. 9). The other set is called the 


Trailing edge 



rotated coordinates and uses a conventional coordinate system for airfoils. In the ro- 
tated coordinate system, the abscissa is tangent to the radii of the leading and trailing 
edges on the pressure side of the blade, and the ordinate is tangent to the leading -edge 
radius. The abscissa is labeled L for length, and the ordinate is labeled H for height 
(see fig. 10). 

The unrotated coordinates for each blade section are calculated by interpolation of 

the blade-element coordinates in the same manner as in the stacking procedure. The 

y and y coordinates, which define the suction and pressure surface profiles of the 
s p 

blade section, are calculated for the complete range of z values. Since the z values 
generally extend beyond both edges of a blade section, a few nonexistent points are cal- 
culated. The leading-edge and trailing-edge coordinates on the suction and pressure 
surfaces of the blade sections are calculated by interpolation of the x-z coordinates of 
the blade elements to obtain the z -coordinates, and then interpolation of the y-z coordi- 
nates of the blade-section surfaces to obtain the y- coordinates. The coordinates of the 
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maximum thickness points and the transition points on the suction and pressure surfaces 
are calculated in the same manner. The center -of -area coordinates are obtained by in- 
terpolation of the x-y and the x-z coordinates of the stacking line. The coordinates of 
the leading -edge, maximum thickness, transition, and trailing -edge points on the center - 
line are obtained by interpolation of the x-y and x-z coordinates of the points on the 
blade elements. 



Figure 10. - Rotated blade section. 


The rotated coordinates of a blade section are calculated by rotation and translation 
of the unrotated coordinates. The angle of rotation y is the angle from the z-axis to the 
L-axis (see fig. 10) and is calculated by equation (53). The rotated coordinates of the 
leading-edge, maximum thickness, transition, and trailing-edge points on the centerline, 
the suction surface, and the pressure surface are directly calculated by equations (54) 
and (55): 


y = sin - * 


^oc ^ic^"V^ Z oc z ic^ + ^oc”^ic^ ( « ' ) - ^ z oc z ic^ 




(z - z. )^ + (y - y. )^ 
v oc ic' vjr oc • y ic' 


(53) 


t. 


H = ( y - y ic )cos y - (z - z. c )sin y + ± 


1C' 


L. 

L = (y - y ic )sin y + (z - z. c )cos y + -i 


1C' 


(54) 

(55) 


The coordinates of the center of area and a reference point, the stacking point of the hub 
blade element, are also calculated by equations (54) and (55). These two points will 
coincide if the stacking line is not tilted. 
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The rotated coordinates of the suction and pressure surface profiles for a blade 
section are obtained at equal increments along the L-axis. These coordinates are cal- 
culated by interpolation of coordinates obtained by rotation and translation of the unro- 
tated coordinates. These coordinates are calculated only for points actually on the 
blade -section surfaces. 

Along with the blade-section rotated coordinates, several parameters which pertain 
to the stress analysis of the blade are calculated. These parameters include the fol- 
lowing: 

(1) Blade -section area, A 

(2) Center -of -area coordinates, L and H 

(3) Moment of inertia about the L-axis, 1^-^ 

(4) Moment of inertia about the H-axis, I HH 

(5) Product of inertia associated with the L-H axes, P HL 

(6) Moment of inertia about L-axis translated to the center of area, Illca 

(7) Moment of inertia about H-axis translated to the center of area, IjjHC A 

(8) Product of inertia associated with the L-H axes translated to the center of 

area, P HLCA 

(9) Angle to the axis of minimum moment of inertia from the L-axis, /3 

(10) Minimum moment of inertia about an axis through the center of area, I m ^ n 

(11) Maximum moment of inertia about an axis through the center of area, I max 
The equations for calculating these parameters are 



- H p )dL 
- H p )dL 



(56) 

(57) 

(58) 

(59) ' 


HH 


L ( H s - H p) dL 


HL 


f L m 

J L=0 

jCTS*"» 


(60) 


(61) 


T LLCA “ *LL “ r2a 


^HCA “ J HH “ l2a 


P HLCA “ P HL " HLA 


(62) 

(63) 

(64) 


(3 = itan" 1 
2 VI 


2P 


HLCA 


HHCA ‘ J LLCA , 


^max *LLCA + WlHCA Wnin 


(65) 


Win ' “ ^LLCA + WmCA^ + \ ^LLCA ~ ^HHCA^ C0S ^ 2 ^^ “ P HLCA sin ( 2 ^) ( 66 ) 


(67) 


The integrals in the preceding equations are evaluated by a numerical integration techni- 
que. 

The form of the output is shown in appendix C for a sample case . 


CONCLUDING REMARKS 

The equations and procedure for defining a complete compressor blade have been 
presented in the preceding sections. Specific details of a computer program which in- 
corporates these equations and procedures are given in appendix C. The details include 
a FORTRAN IV source deck listing of the program, definitions of the program variables, 
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descriptions of the subroutines, the input format, and an output listing for a sample 
blade. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, June 26, 1969, 
720-03-00-64-22. 
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APPENDIX A 
SYMBOLS 


A 

C 

f (k,K 1 ,Z,R,R 1 ) 

H 

I 

L 

R 

r 

s 

d 

t 

X 


y 

z 

a 

0 

r 

6 

€ 

C 

5? 


K 


area of blade section 
rate of turning, -dfc/ds 

function describing relation of e - to k and R for given values 
of Rp Kj, and £ 

height coordinate for blade section 

blade section moment of inertia 

length coordinate for blade section 

distance from vertex to point on cone 

radial coordinate in cylindrical system 

path length along blade-element centerline or surface 

path length along blade thickness line 

blade thickness 

distance from axis of rotation along radial line passing through hub- 
element stacking point (fig. 8) 

coordinate perpendicular to z in constant x-plane (fig. 8) 

axial coordinate from hub-element, leading-edge center 

cone half- angle (fig. 1) 

angle of axis of minimum moment of inertia to L-axis (fig. 10) 

angle of axis of rotation to L-axis (fig. 10) 

circumferential angle coordinate of stacking line 

angular coordinate on conic surface as measured from ray passing 
through blade -element, leading-edge center (fig. 1) 

convenient constant on a segment, eq. (7) 

lean angle of stacking line in r-9 plane (fig. 8) (positive in positive 
9- direction) 

local blade angle, the angle between the local R and the tangent to 
the local blade-element centerline or surface path (fig. 1) 
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X 


lean angle of stacking line in r-z plane (fig. 8) (positive in positive z- 
direction) 

9 circumferential angle coordinate in cylindrical coordinate system (positive 9 - 

direction is from pressure surface to suction surface) 

Subscripts: 

c blade centerline 

ca center of area 

h hub element 

i inlet segment or leading edge 

j index denoting axial location 

m maximum thickness point 

max maximum value 

min minimum value 

n normal to blade-element centerline 

0 outlet segment or trailing edge 

p pressure surface 

s suction surface 

sp stacking point 

t transition point between segments of blade 

1 arbitrary reference or known value 

2 known value 

3 known value 

Superscript: 

center-of-area coordinate 
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APPENDIX B 


PARTICULAR FORMS OF THE GENERAL EQUATIONS 
Solutions of R, R^) 

For a line of constant turning rate C on a conical surface, the differentials of the 
radial and angular coordinates can be expressed as follows: 

dR = _cosjc dK (Bl) 

C 


and 


de = - 


sin k 
RC 


dK 


(B2) 


Equation (Bl) integrates to 


R - Rh = — (sin k, - sin k) (B3) 

1 C 1 

Rearrangement of equation (B3) yields a characteristic constant £ for a line of con- 
stant C 


£ = RC + sin k = RjC + sin k j (B4) 

By solving equation (B4) for RC and replacing RC in equation (B3), de becomes a 
function of k 


de = sin - - dK (B5) 

sin k - ? 

However, if C = 0, k is a constant, and equation (B5) is indeterminate. A different 
equation is required for this special case. For C = 0 or for constant k, e is a func- 
tion of R, and the differential can be expressed as follows: 

de = tan /c 1 — (B6) 

1 R 
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In general, the indefinite integral of de is given by 


e - € 1 =f (k, 


(B7) 


where the function f(/c, k j, ?,R,R^) has four different solutions dependent on k, 
and The forms of the function are as follows: 

(1) If k — (i.e., C = 0), 


i(n, k^, C, R, Rj) = tan k In 



(B8) 


(2) If k * /Cj and > 1, 


Hk,K±, ?,R, Rj) = K - fCj + 







r _ 


” hr \ ~ 


-4 

1 “ 5 tan fl) 

H 

1 ' ? 


^tan” 

V 2 / 

- tan " 1 

\2 / 

V 

V? 2 - 1 

S 2 - 1 


V 



*” “ 

J 


(B9) 


(3) If k # k x and £ 2 < 1, 

f^j/Cjj^R, R x ) = K - k 1 


- 

1 - t; tan ^ - £ 2 


1 - 5 tan - Vl - ? 2 

o * 


In 

2 

- Zn 

A 


1 - £ tan £ + Vl - 

9 T 

1 - ^ tan — + - £ 2 

2 



Ci 



- 


(4) If k * Kj and £ = ±1, 

f Oe, k v ^,R,R 1 )=k-k 1 ± [tan^- ± ^ - tan^- ± ^ 


(BIO) 


(Bll) 
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Equations for Inlet Segment of Centerline 

The equations for the inlet segment of the centerline are derived from equations 
(B4) and (B7) with the appropriate constants 


he = R ic C ic + sin *ic 

(B12) 

K c = sm ^ic " C ic R c^ 

(B13) 

e c - e ic + f ^c’ ^ic’ ^ic’ R c> R ic^ 

(B14) 


These equations apply for s R^. For convenience, the center of the leading edge is 
used as a reference and, thus, e^ c = 0. 

Equations for Outlet Segment of Centerline 

The equations for the outlet segment of the centerline have the same form as those 
for the inlet segment centerline, but have different C and £ constants 


?oc = R tc C oc + sin K tc 

(B15) 

K c = sin_1 ^oc - C oc R c) 

(B16) 

e c = e tc +f ( K c’ K tc’W R c’ R tc ) 

(B17) 

where e <;c is evaluated at the end of the inlet segment centerline or 
point by equation (B14). These equations apply for R c > R^. 

centerline transition 

Surface Coordinates at Ends of Thickness Path 


In the R-e coordinate system, a thickness path is described by a line of constant 
angle k , which is perpendicular to the centerline 


(B18) 
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In equation (B18), the angle to the suction surface is given by the plus sign, and the 
angle to the pressure surface is given by the minus sign. 

The differential equations for the thickness path in terms of the path direction k 
and the path distance a are 


dR = cos /c n da 


(B19) 


and 


Rde = sin « n da 


(B20) 


Integration of equation (B19) gives 


R - R = a cos k 
c n. 

Substitution of da from equation (B19) into equation (B20) gives 

de = tan /c n — 
n R 


(B21) 


(B22) 


Integration of equation (B22) gives 


e 


= tan 



(B23) 


However, if k q = 0 , K n = ±7r/2 and R = R c , and equation (B23) becomes indeterminate. 
Since R is a constant for this special case, equation (B20) can be integrated as 


e 


sin k 


e c = 


n 


R, 


(B24) 


where k = ±jt/2. 

Blade thickness is specified at three locations: the leading edge, the maximum 
thickness point, and the trailing edge. At these three locations, the suction surface and 
pressure surface coordinates are calculated by the use of the appropriate thickness 
value and corresponding blade centerline angle in equations (B21) and either (B23) or 
(B24). 
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On the suction surface, 


r -» (i, m o) . 

R (i,m,o)s R (i,m,o)c g sin K c 


"(i,m,o)s fc (i,m,o)c 


+ cot K( , In 

R (i,m,o)s 


e (i,m,o)p e (i,m,o)c + 


2R (i,m,o)c 


On the pressure surface, 


t/* \ — R/ • \ - — - — sin k 

(i,m,o)p (i,m,o)c 2 c 


€. / • \ — ^ / • \ + cot k In 

fc (i,m,o)p (i,m,o)c c 


R (i,m,o)c 

R (i,m,o)p 


Or, if k c = 0, 


= e _ *(*. ”*.<>) 

(i,m,o)s (i,m,o)c 2R 

(i, rn, o)c 
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APPENDIX C 


DESCRIPTION OF COMPUTER PROGRAM 

The blade coordinate computer program incorporates the equations and calculation 
procedures presented in this report to compute the cross-section coordinates of a com- 
pressor blade composed of multiple -circular -arc elements on conical surfaces. In ad- 
dition to the coordinates, parameters for stress analysis (such as area, center of area, 
and moments of inertia) are also computed. The program consists of a main program 
and several subprograms. It is written in FORTRAN IV. The run time on a direct- 
coupled IBM 7044-7094 system is approximately 0. 01 minute per given blade element. 

The information in the following sections is intended to aid in the use of the program 
and in the understanding of its logic. Included are a description of the input, definitions 
of program variables, descriptions of subprograms, a listing of the program, and a 
sample output. 


Description of Input 

The format for the input cards is shown in table I. The first card in a set of data 
is the title card. It is used to identify the data with alphanumeric information, which is 
printed out with the output data. The second card is a general card for specification of 
single -value variables. The definitions of these variables are as follows: 

ETA tangential lean angle of stacking line rj, in degrees (positive in direction 

from pressure surface toward suction surface) 

LAMDA axial lean angle of stacking line X, in degrees (positive in direction from 
inlet toward outlet) 

XNR number of blade elements 

OP1 n um ber of specified radial locations for desired blade sections (If none are 

specified (i.e. , OP1 = 0. 0), program computes blade sections at radial 
locations of stacking points for all blade elements. ) 

OP2 control variable for output of calculated blade-element parameters (angles 

and turning rates) and coordinates (Blade -element output is printed out 
if OP2 = 1. 0.) 

TNLMT tolerance limit for blade-element stacking iteration (If the tolerance limit 
is set too small, the stacking procedure will require an excessive 
number of iterations and may not converge. ) 
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The next set of cards specifies the geometry of the blade elements. As shown in table I 
for the first variable, RI, data for each variable begins in the first data space on a card 
and continues in succeeding spaces and cards for a total of XNR spaces . The maximum 
number of data per variable (i.e. , number of blade elements) is 24. The definitions of 
the input blade-element variables are given in the main text under the section DEVELOP- 
MENT OF EQUATIONS FOR BLADE-ELEMENT LAYOUT. The correspondence between 
variable names and variable symbols is as follows: 

RI inlet radius, r. 

lv 

RO outlet radius, r 

Uv 

TI inlet blade thickness, t. 

TM maximum blade thickness, t m 

TO outlet blade thickness, t 

KIC blade centerline angle at inlet, K ic 

KTC blade centerline angle at transition point, k^ c 

KOC blade centerline angle at outlet, k qc 

ZMC axial distance to maximum thickness point from inlet, z mc - z^ c 

ZTC axial distance to transition point from inlet, z^ c - z. £ 

ZOC axial distance to outlet from inlet, z QC - z^ c 

The last set of cards specifies the radial locations XQ of the desired blade sections. 

The XQ input begins in the first data space and continues in succeeding spaces for a total 

of OP1 spaces. The maximum number of blade sections is 24. 

All input data are floating-point numbers. All input angles are in degrees. All 

variables with length dimensions must have the same unit length. The inlet and outlet 

radii r. and r of a blade element cannot be identical. A difference in radii of at 
X c oc 

least 0. 1 percent of the axial length of the blade element is recommended. 


Main Program Variables and Definitions 

The following is a list of program variable names with their corresponding symbols 


or definitions. 



FORTRAN 

Mathematical 

Definition 

variable 

symbol 


AREA 

A 

Area of blade section 
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FORTRAN 

Mathematical 

Definition 

variable 

symbol 


ALP(I) 

O' 

Cone half-angle 

BETA 

/3 

Angle of I min axis (eq. (65)) 

BIC(I) 

^1C 

Blade centerline, inlet segment constant (eq. (7) 

BIP(I) 

s ip 

Pressure surface, inlet segment constant (eq. (7)) 

BIS (I) 

S 1S 

Suction surface, inlet segment constant (eq. (7)) 

BOC(I) 

b oc 

Blade centerline, outlet segment constant (eq. (7)) 

BOP(I) 

b op 

Pressure surface, outlet segment constant (eq. (7)) 

BOS(I) 

e 

b os 

Suction surface, outlet segment constant (eq. (7)) 

CAPPA 

K 

Local blade angle (eq. (23)) 

CAPRZ 

R 

Distance from vertex to point on cone (eq. (14)) 

CAS (I) 

C os 

Suction surface, outlet segment, rate of turning 
(eq. (15)) 

CIC(I) 

C ic 

Blade centerline, inlet segment, rate of turning 

CIP(I) 

c ip 

Pressure surface, inlet segment, rate of turning 

CIS (I) 

C is 

Suction surface, inlet segment, rate of turning 

COC(I) 

C oc 

Blade centerline, outlet segment, rate of turning 

COP(I) 

C op 

Pressure surface, outlet segment, rate of turning 

CRCG(I) 

R 

sp 

Stacking point radius 

DZ 

Az 

Axial coordinate increment 

EIC 

6 ic 

Blade centerline, inlet segment, angular coordi- 
nate (EIC = 0.0) 

EIP(I) 

€ iP 

Pressure surface, inlet segment, angular coordi- 
nate (eq. (11)) 

EIS(I) 

e is 

Suction surface, inlet segment, angular coordinate 

EMC (I) 

e mc 

Blade centerline, maximum thickness point, an- 
gular coordinate 

EMP(I) 

e mp 

Pressure surface, maximum thickness point, an- 
gular coordinate 
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FORTRAN 

Mathematical 

Definition 

variable 

symbol 


EMS(I) 

e ms 

Suction surface, maximum thickness point, angular 
coordinate 

EOC(I) 

e oc 

Blade centerline, outlet segment, angular coordi- 
nate 

EOP(I) 

e op 

Pressure surface, outlet segment, angular coordi^ 
nate 

EOS (I) 

6 os 

Suction surface, outlet segment, angular coordinate 

EPC 

e c 

Blade centerline, angular coordinate on z cuts 

EPP 

e p 

Pressure surface, angular coordinate on z cuts 

EPS 

e s 

Suction surface, angular coordinate on z cuts 

ETA 

V 

Lean angle of stacking line in r-0 plane 

ETC(I) 

e tc 

Blade centerline, transition point, angular coordi- 
nate 

ETP(I) 

e tp 

Pressure surface, transition point, angular coor- 
dinate 

ETS(I) 

e ts 

Suction surface, transition point, angular coordi- 
nate 

GAMX(K) 

y 

Angle of L-axis from axis of rotation (eq. (53)) 

HBAR 

H 

Blade section, center-of-area coordinate (eq. (58)) 

I 

— 

Index used to denote blade element 

ICASE(I) 

— 

Integers 1, 2, or 3 denoting whether transition is 
ahead of, equal to, or behind maximum thickness 

IHH 


Blade section moment of inertia (eq. (60)) 

IHHCG 

^HCA 

Blade section moment of inertia (eq. (63)) 

ILL 

J LL 

Blade section moment of inertia (eq. (59)) 

ILLCG 

*LLCA 

Blade section moment of inertia (eq. (62)) 

IMAX 

I 

max 

Blade section, maximum moment of inertia 
(eq. (67)) 

IMIN 

■'min 

Blade section, minimum moment of inertia 
(eq. (66)) 
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FORTRAN 

variable 


Mathematical 

symbol 


Definition 


J 

j 

K 

— 

KIC(I) 

K ic 

KIP(I) 

K iP 

KIS(I) 

“is 

KM (I) 

K m 

KOC(I) 

*oc 

KOP(I) 

“op 

KOS(I) 

K os 

KTC(I) 

K tc 

KTP(I) 

K tp 

KTS(I) 

“ts 

LAMDA 

X 

NR 

— 

NXQ 

— 

NZ 

— 

OP1 

— 

OP2 

* 

— 

PHL 

P HL 

PHLCG 

P HLCA 

RCG(I) 

R sp 

RI(I) 

r ic 

RIC(I) 

R ic 

RIP(I) 

R ip 

RIS(I) 

R is 


Index denoting position in z-direction 

Index denoting position in x-direction 

Blade centerline, leading-edge, local blade angle 

Pressure surface, leading-edge, local blade angle 

Suction surface, leading-edge, local blade angle 

Maximum thickness point, local blade angle 

Blade centerline, trailing-edge, local blade angle 

Pressure surface, trailing-edge, local blade angle 

Suction surface, trailing-edge, local blade angle 

Blade centerline, transition point, local blade 
angle 

Pressure surface, transition point, local blade 
angle 

Suction surface, transition point, local blade angle 
Lean angle of stacking line in r-z plane 
Number of input radii 
Number of blade sections 
Number of z values 

Number of blade-section locations specified in input 

Control variable for printed output of blade- 
element coordinates and parameters 

Blade section product of inertia (eq. (61)) 

Blade section product of inertia (eq. (64)) 

Stacking point radius 

Blade centerline, leading-edge, radial coordinate 
Blade centerline, leading-edge radius 
Pressure surface, leading-edge radius 
Suction surface, leading-edge radius 
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FORTRAN 

variable 


Mathematical 

symbol 


Definition 


RMC(I) 


RMP(I) 

R mp 

RMS (I) 

R ms 

RO(I) 

r oc 

ROC (I) 

R„„ 

oc 

ROP(I) 

R op 

ROS(I) 

R os 

RTC(I) 

R tc 

RTP(I) 

R tp 

RTS (I) 

R ts 

T,T1,T2,T3, 

— 

T4, T5 


THECG(I) 

0 sp " 

THETA (I) 

e ic 

THETAC(I, J) 

6 c 

THE TAP (I, J) 

0 P 

THETAS(I, J) 

0 s 

TI(I) 


TM(I) 

fc m 

TNLMT 

— 

TNORM1 

S 

TO(I) 

*0 

Vi,V 2 ,V 3 ,V 4 

— 

X(K) 

X 

XCG(I) 

X„„ 

sp 

XHCG 

— 

XHIC 




Blade centerline, maximum thickness point radius 
Pressure surface, maximum thickness point radius 
Suction surface, maximum thickness point radius 
Blade centerline, trailing-edge, radial coordinate 
Blade centerline, trailing-edge radius 
Pressure surface, trailing-edge radius 
Suction surface, trailing-edge radius 
Blade centerline, transition point radius 
Pressure surface, transition point radius 
Suction surface, transition point radius 
Temporary storage locations 

Relative stacking point, circumferential angle co- 
ordinate 

Blade centerline, leading-edge, circumferential 
angle coordinate 

Blade centerline, circumferential angle coordinate 
Pressure surface, circumferential angle coordinate 
Suction surface, circumferential angle coordinate 
Leading-edge blade thickness 
Blade thickness at maximum thickness point 
Blade-element stacking tolerance limit 
Blade-element stacking tolerance (eq. (99)) 
Trailing-edge blade thickness 
Temporary storage locations 

Computed values of x-coordinate for blade sections 
Stacking point x- coordinates 
Blade section, center-of-area H-coordinate 
Blade section, centerline, leading-edge H-coordinate 
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F 


FORTRAN 

variable 

Mathematical 

symbol 

Definition 

XHIP 

— 

Blade section, pressure surface, leading-edge H- 
coordinate 

XHIS 

— 

Blade section, suction surface, leading-edge H- 
coordinate 

XHMC 

— 

Blade section, centerline, maximum thickness 
point H- coordinate 

XHMP 

— 

Blade section, pressure surface, maximum thick- 
ness point H- coordinate 

XHMS 

— 

Blade section, suction surface, maximum thickness 
point H- coordinate 

XHOC 

H oc 

Blade section, centerline, trailing-edge H- 
coordinate 

XHOP 

H op 

Blade section, pressure surface, trailing-edge H- 
coordinate 

XHOS 

H os 

Blade section, suction surface, trailing-edge H- 
coordinate 

XHTC 

H tc 

Blade section, centerline, transition point H- 
coordinate 

XHTP 

H tp 

Blade section, pressure surface, transition point 
H- coordinate 

XHTS 

H ts 

Blade section, suction surface, transition point FI- 
coordinate 

XIC(I) 

x ic 

Blade centerline, leading-edge x-coordinate 

XIP(I) 

x ip 

Pressure surface, leading-edge x-coordinate 

XIS(I) 

x is 

Suction surface, leading-edge x-coordinate 

XLCG 

— 

Blade section, center-of-area L-coordinate 

XLIC 

L ic 

Blade section, centerline, leading- edge L- 
coordinate 

XLIP 

L ip 

Blade section, pressure surface, leading-edge L- 
coordinate 


I 
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FORTRAN 

variable 


Mathematical 

symbol 


Definition 


XLIS 

L is 

Blade section, suction surface, leading-edge L- 
coordinate 

XLMC 

L mc 

Blade section, centerline, maximum thickness 
point L- coordinate 

XLMP 

^mp 

Blade section, pressure surface, maximum thick- 
ness point L- coordinate 

XLMS 

L ms 

Blade section, suction surface, maximum thickness 
point L- coordinate 

XLOC 

L oc 

Blade section, centerline, trailing-edge L- 
coordinate 

XL OP 

L op 

Blade section, pressure surface, trailing-edge L- 
coordinate 

XLOS 

L os 

Blade section, suction surface, trailing-edge L- 
coordinate 

XLSP 

— 

Blade section, reference point (hub blade element 
stacking point) L- coordinate 

XLTC 

L tc 

Blade section, centerline, transition point L- 
coordinate 

XLTP 

L tp 

Blade section, pressure surface, transition point 
L- coordinate 

XLTS 

L ts 

Blade section, suction surface, transition point 
L- coordinate 

XMC(I) 

x mc 

Blade centerline, maximum thickness point x- 
coordinate 

XMP(I) 

x mp 

Pressure surface, maximum thickness point x- 
coordinate 

XMS(I) 

x ms 

Suction surface , maximum thickness point x- 
coordinate 

XOC(I) 

x oc 

Blade centerline, trailing-edge x- coordinate 

XOP(I) 

x op 

Pressure surface, trailing-edge x-coordinate 

XOS(I) 

x os 

Suction surface, trailing-edge x-coordinate 
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FORTRAN 

variable 


Mathematical 

symbol 


Definition 


XQ(K) 

— 

XP(I,J) 

— 

XTC(I) 

x tc 

XTP(I) 

X tp 

XTS(I) 

x ts 

YIC(I) 

yic 

YIP(I) 

y*p 

YIS(I) 

^is 

YMC (I) 

Vc 

YMP(I) 

y mP 

YMS(I) 

^ms 

YOC(I) 

y oc 

YOP(I) 

y 0P 

YOS(I) 

y os 

YP(I,J) 

y P 

YS(I,J) 

y s 

YTC(I) 

ytc 

YTP(I) 

y tp 

YTS(I) 

yts 

YICX(K) 

— 

YIPX(K) 

— 

YISX(K) 

— 

YMCX(K) 

— 

YMPX(K) 

— 


Input values of x-coordinate for desired blade 
sections 

Blade pressure surface abscissa 

Blade centerline, transition point x- coordinates 

Pressure surface, transition point x- coordinates 

Suction surface, transition point x- coordinates 

Blade centerline, leading-edge y-coordinates 

Pressure surface, leading-edge y-coordinates 

Suction surface, leading-edge y-coordinates 

Blade centerline, maximum thickness point y- 
coordinates 

Pressure surface, maximum thickness point y- 
coordinates 

Suction surface, maximum thickness point y- 
coordinates 

Blade centerline, trailing-edge y-coordinates 
Pressure surface, trailing-edge y-coordinates 
Suction surface, trailing-edge y-coordinates 
Blade pressure surface y-coordinates 
Blade suction surface y-coordinates 
Blade centerline, transition point y-coordinates 
Pressure surface, transition point y-coordinates 
Suction surface, transition point y-coordinates 
Value of y^ c at a given blade section 
Value of y. at a given blade section 
Value of y,_ at a given blade section 
Value of y mc at a given blade section 
Value of y mp at a given blade section 


I 


43 



FORTRAN 

variable 


Mathematical 

symbol 


Definition 


YMSX(K) 

— 

Value of y at a given blade section 

ms 

YOCX(K) 

— 

Value of y at a given blade section 

o c 

YOPX(K) 

— 

Value of y Q p at a given blade section 

YOSX(K) 

— 

Value of y Qg at a given blade section 

YTCX(K) 

— 

Value of yj c at a given blade section 

YTPX(K) 

— 

Value of y^ at a given blade section 

YTSX(K) 

— 

Value of y^ s at a given blade section 

ZCG(I) 

z sp 

Stacking point axial coordinate 

ZIC(I) 

z ic 

Blade centerline, leading-edge z-coordinates 

ZIP(I) 

Z ip 

Pressure surface, leading-edge z-coordinates 

ZIS(I) 

z is 

Suction surface, leading-edge z-coordinates 

ZMC(I) 

z mc 

Blade centerline, maximum thickness point z- 


coordinates 

ZMP(I) 

z mp 

Pressure surface, maximum thickness point z- 


coordinates 

ZMS(I) 

z ms 

Suction surface, maximum thickness point z- 


coordinates 

ZOC(I) 

z oc 

Blade centerline, trailing-edge z-coordinates 

ZOP(I) 

z op 

Pressure surface, trailing-edge z-coordinates 

ZOS(I) 

z os 

Suction surface, trailing-edge z-coordinates 

ZTC(I) 

z tc 

Blade centerline, transition point z-coordinates 

ZTP(I) 

z tp 

Pressure surface, transition point z-coordinates 

ZTS(I) 

z ts 

Suction surface, transition point z-coordinates 

ZX(J) 

z. 

J 

Values of equally spaced z- increments computed 


to obtain x-y cuts 

ZICX(K) 

— 

Value of z^ c at a given blade section 

ZIPX(K) 

— 

Value of z^ at a given blade section 

ZISX(K) 

— 

Value of z. at a given blade section 

IS 
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FORTRAN 

Mathematical 



Definition 

variable 

symbol 




ZMCX(K) 

— 

Value of 

z mc 

at a given blade section 

ZMPX(K) 

— 

Value of 

z mp 

at a given blade section 

ZMSX(K) 

— 

Value of 

z ms 

at a given blade section 

ZOCX(K) 

— 

Value of 

z oc 

at a given blade section 

ZOPX(K) 

— 

Value of 

Z op 

at a given blade section 

ZOSX(K) 

— 

Value of 

z os 

at a given blade section 

ZTCX(K) 

— 

Value of 

z tc 

at a given blade section 

ZTPX(K) 

— 

Value of 

z tp 

at a given blade section 

ZTSX(K) 

— 

Value of 

z ts 

at a given blade section 


Description of Subroutines 

The subroutines used in this program are listed below along with their call se- 
quence, purpose, and variable definitions. 

Subroutine ITER(K2, C, B, Kl, El, Rl, E2, R2, XK). - A routine to iteratively solve for 
the equation of a constant d/c/ ds curve which passes through two known points and at a 
given slope at one of the points. Refer to equations (7), (11), and (15) for the functional 
relations. 


K2 

*2 

Unknown slope at point 2 

C 

C 

Unknown curvature constant 

B 

K 

Unknown curve constant 

Kl 

K 1 

Known slope at point 1 

El 

e l 

Angular coordinate of point 1 

Rl 

Rl 

Radial coordinate of point 1 

E2 

e 2 

Angular coordinate of point 2 

R2 

*2 

Radial coordinate of point 2 

XK 

— 

An initial estimate of ^ 
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Subroutine ITER1(KT, RT, ET, KM, RM, EM, B, C, RTC, ETC, KTC) . - A routine to 
iteratively solve for the R-e coordinates of the transition point on either the pressure 
surface or the suction surface. Refer to equations (11), (21) to (23), and (25) for the 
functional relations. 


FORTRAN 

variable 

Mathematical 

symbol 

Definition 

KT 

K t(p,s) 

Unknown surface transition point blade angle 

RT 

^(P, s) 

Unknown surface transition point radial coordinate 

ET 

€ t(p, s) 

Unknown surface transition point angular coordi- 
nate 

KM 

K m(p, s) 

Surface maximum thickness point blade angle 

RM 

R m(p, s) 

Surface maximum thickness point radial coordi- 
nate 

EM 

e m(p, s) 

Surface maximum thickness point angular coordi- 
nate 

B 

^(i, o)(p, s) 

Surface curve constant 

C 

C (i,o)(p, s) 

Surface curvature constant 

RTC 

R tc 

Blade centerline, transition point, radial coordinate 

ETC 

e tc 

Blade centerline, transition point, angular coordi- 
nate 

KTC 

K tc 

Blade centerline, transition point, blade angle 

Subroutine SINTP(Z, W, N, XI, Y2). - A routine which uses a second-order Lagrangian 

algorithm for interpolation and a linear extrapolation method. Refer to equations (40) 
to (43) for the functional relations. 

Z 

X j(p,s) 

Abscissa vector 

w 

y j(p, s) 

Corresponding ordinate vector 

N 

— 

Number of points in the given vector 

XI 

X 

Given argument 

Y1 

y (p, s) 

Interpolated ordinate, (i. e. , Y1 = W(X1)) 
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Subroutine CGS(YCG, ZCG, YIPX, ZIPX, YISX, ZISX, YPX, YSX, NZ, ZX, YOPX, ZOPX, 
YOSX, ZOSX). - A routine to calculate the center of area of a blade section. 


FORTRAN 

variable 

Mathematical 

symbol 

Definition 

YCG 

— 

Value of y of the blade section 

ca 

ZCG 

— 

Value of z„„ of the blade section 
ca 

YIPX 

— 

Value of y. of the blade section 

ZIPX 

— 

Value of of the blade section 

YISX 

— 

Value of y. _ of the blade section 

ZISX 

— 

Value of z. of the blade section 

IS 

YPX(J) 

y P 

Blade pressure surface ordinates of the blade 
section 

YSX(J) 

y s 

Blade suction surface ordinates of the blade sec- 
tion 

NZ 

— 

Number of z stations 

ZX(J) 

z. 

J 

Values of z for all z stations 

YOPX 


Value of y Q p of the blade section 

ZOPX 

— 

Value of z Q p of the blade section 

YOSX 

— 

Value of y_ of the blade section 

ZOSX 

— 

Value of z of the blade section 

OS 

Subroutine FIX(A, B). - A routine to determine the arcsin of a given value and pre- 
vent computation of the arcsin of a value greater than 1. 0 or less than -1. 0. 

A 

— 

Given value 

B 

— 

Computed arcsin (A) 

Subroutine CGS1(X). - A routine to calculate the center of area, moment of inertia, 
minimum moment of inertia, and axis of minimum moment of inertia for a blade section 
in the L-H coordinate system. 

X 

L 

Chordwise abscissa 

HP 

H P 

Blade pressure surface ordinate 

HS 

H s 

Blade suction surface ordinate 
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FORTRAN 

variable 

Mathematical 

symbol 

Definition 

N 

— 

Number of chordwise abscissas 

XI 

— 

Chordwise abscissa 

V 

— 

Blade suction surface ordinate associated with XI 

VI 

— 

Blade pressure surface ordinate associated with 
XI 

V2 

— 

Temporary storage for function value 

V3 

— 

Temporary storage for integral value 

AREA 

A 

Blade -section area 

LBAR 

L 

Center -of -area coordinate 

HBAR 

H 

Center-of-area coordinate 

IHH 

I HH 

Moment of inertia about H-axis 

ILL 

*LL 

Moment of inertia about L-axis 

PHL 

P HL 

Product of inertia, P HL = jy HL dA 

Subroutine 

XMAX(X, XM,N). 

- A routine which selects the maximum value of a vec- 

tor. 



X 

— 

Given vector 

XM 

— 

Maximum value of X 

N 

Subroutine XMIN(X,XM,N). - 

tor. 

Number of points in X 

A routine which selects the minimum value of a vec- 

X 

— 

Given vector 

XM 

— 

Minimum value of X 

N 

— 

Number of points in X 

Subroutine NEED(IC,I). - A routine which directs the computation of the curves for 
the pressure and suction surfaces of a blade element. 

IC 


Integer (1, 2, or 3) denoting whether transition 
point is at, ahead of, or behind maximum thick- 
ness point 

I 

— 

Index corresponding to blade element to be com- 


puted 
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Function SUBF(X, XO, B, R, RO). - A subprogram to compute the function 
f (k, k q , £, R, R q ) as given in appendix B. 


FORTRAN 

Mathematical 

Definition 

variable 

symbol 


X 

fC 

Slope of curve at a point 

XO 

K o 

Slope of curve at a reference point 

B 

? 

Curve constant 

R 

R 

Radial coordinate of point 

RO 

R o 

Radial coordinate of reference point 


Subroutine INTGR(L, X1,X2,X3). - A routine to evaluate three definite integrals de- 
fined by equations (28), (29), and (30). 


L 


— 

Index denoting blade element 

XI 


— 

Value of integral J' e dA 

X2 


— 

Value of integral J' R dA 

X3 


— 

Value of integral f dA 


Function ADJ(D) 

. - A subprogram to adjust the increment D to a value not less 

than 

D and having a 

single significant figure of 1, 2, or 5. 

D 


— 

Increment 


Subroutine RAEP(RP, RS, EP, ES, RC, EC, XKC, TC). - A routine to calculate the 

conical coordinates of surface points at the ends of a thickness path of a blade element 

Refer to equations (16), (19), (21), 

and (22) for functional relations. 

RP 


R (i,m,o)p 

R-coordinate for pressure surface point 

RS 


R (i,m,o)s 

R-coordinate for suction surface point 

EP 


€ (i, m , o)p 

e-coordinate for pressure surface point 

ES 


€ (i, m , o)s 

e-coordinate for suction surface point 

RC 


R (i, m, o)c 

R-coordinate for centerline point 

EC 


e (i,m,o)c 

e-coordinate for centerline point 

XKC 

K (i,m,o)c 

k at centerline point 

TC 


t (i,m,o) 

Blade thickness 
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Subroutine FNTGRL(N,DX, FX, SFX). - A Lewis system subroutine to numerically 
evaluate the integral of a function defined at any number of equally spaced intervals. 

FORTRAN 

variable 

Mathematical 

symbol 

Definition 

N 

— 

Number of stations 

DX 

dx 

Size of interval 

FX 

f(x) 

Values of function at each station 

SFX 

f(x)dx 

Jo 

Values of integral 

Subroutine 
values in the X- 
to maintain the 

SORTXY(X, Y,N). - A Lewis system subroutine to rearrange the N 
-array in order of increasing size and move the values of the Y-array 
original pair relations. 

X 

— 

Independent array 

Y 

— 

Dependent array 

N 



Number of values 


FORTRAN IV Source Deck Listing 


C BLADE COORDINATE PROGRAM FOR BLADE WITH 

C 2 PARTS OF DIFFERENT CURVATURE 

C 

C OP1 ENTER OWN X COORDS 

C OP2 X-Y-Z CONICAL SECTION COORDINATES 

DIMENSION ALP(24), BIC(24), BIP(24>, BIS(24), BOC ( 24 ) , B0P(24), BO 
1S(24), CIC(24), Cl P( 24) » CTS(24), COCI24), C0PI24), CAS(24), DELTI 
224), EIPI24), EISI24), EMCI24), EMPJ24), EMSI24), F0C.(24), E0PI24) 
3, E0SC24), ETC! 24) , ETP(24>, ETSI24), ICASEI24), KIC(24), KIP(24>, 
4 K I S ( 24 ) » KM (24 ) , KOC(24), K0P(24), KOS ( 24 ) , KTC(24), KTP(24), KTS 
5(24), RCG( 24 ), RI(24), RIC(24), RIP(24>, RIS(24), RMC(24), RMP(24) 
6, RMS( 24) , RCC ( 24) , R0P(24), R0S(24), RTC(24), RTP(24), RTS(24), T 
7HECG124), THETA( 24) , T I ( 24) , TM(24), T0<24), XCG(24), XIC(24), XIP 
8(24), XIS( 24 ), XMC ( 24) , XMP(24), XMSC24), XOC(24), XOP(24), XOS(24 
9), XTC ( 24 ) , XTP ( 24 ) , XTS(24), YCG(24), YIC(24), YIP(24), YIS(24), 
$Y MC ( 24 ) , YMP(24), YMSI24), Y0C(24), Y0P(24), Y0S(24), YTC(24), YTP 
$(24), YTS ( 24 ), ZCG ( 24 ) , ZIC(24), ZIP(24), ZIS(24), ZMC(24), ZMP(24 
$) , ZMS ( 24 ) , ZOC ( 24 ) , Z0P(24>, ZOSJ24), ZTC(24), ZTP(24), ZTS(24) 
DIMENSION Z X (3 2 ) , YPXI32), YSX(32), XP(32,24), XS(32,24), YP(32,24 
1), YS( 32,24), THETAC ( 32,24) , THETA P ( 32 ,24) , THETAS < 32,24) 

DIMENSION V ( 56 ) , Vl(56), V2<56), V3(56), V4(56), V5(56) 

DIMENSION X ( 24 ) , XQ(24), GAMX(24), TIX(24), TMX(24), T0X(24), YCGX 
1(24), YICX( 24) , YI PX( 24) , YISX(24), YMCXI24) , YMPX(24), YMSXI24), 
2Y OCX ( 24 ) , YOPX ( 24) , Y0SX(24), YTCX<24), YTPX(24), YTSXI24), ZCGX(2 
34), ZICXI24), Z IPX ( 24) , ZISX(24), ZMCXI24) , ZMPXI24), ZMSX ( 24 ) , ZO 
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c 

1 


2 

C 


3 


4CX( 2.4) » ZOPX 12 4 ) » ZOSX ( 24 ) • ZTCX(24) , ZTPXC24), ZTSX(24) 

DIMENSION XL (56 ) , XHP(56), XHSI56) 

DIMENSION RO ( 24 ) » ZI<24). CRCG(24), TITLEI12) 

DIMENSION TALPI24), SALP(24), CALP (24) 

COMMON ALP.BOC, BIC»B IP.BI S,BOP ,BOS ,CI P ,C I S ,COP ,CAS ,C IC , El P , El S , EMC 
U EMP.EMS, EOC ,EOP ,EOS , ETC , ETP ,E TS, KI C , KI P ,KI S ,KM, KOC , KOP , KOS , KTC , KT 
2P»KT£»RI,RO»RIC ,RIP ,RI S.RMC ,RM P,RMS, ROC , ROP ,ROS , RTC, RTP, RTS ,T I ,TM, 
3T 0, ZNC, ZOC » Z TC » NR, COC 

REAL KIC* KIPyKI S, KM, KOC, KOP » KOS, KTC ,KTP,KTS, LAMDA 
REAL LBARylLLylHHyl LLCGy IHHCGy IMINylMAX 

COMMON /EXTRA/ XL y XHSy XHP y LOUT y LBAR y HBAR y I LL • I HH , PHL , ARE A t I LLCGy IH 
lHCGyPHLCG, 8ETA,IMIN» IMAX 
DATA. DEGRAD/57.2958/ 

INPUT 


READ 15,43) TITLE 

READ 15,44) ET A, LAMDA, XNR ,0P1 ,0P2, TNLMT 

REWLND 2 

NR=XNR 

R FAD (5,44) <RI(I)yI=lyNR) 

REAO (5,44) (R0( I ) ,I=1»NR) 

READ (5,44) (TI ( I ) y 1 = 1 ,NR ) 

READ (5,44) (TM(I),I=1,NR) 

R FAD (5,44) (T0(I),I=1,NR) 

READ (5,44) (K IC ( I ) , I =1 , NR ) 

READ (5,44) (K TC ( I ) , 1 = 1 , NR ) 

READ (5,44) (KOC (I ) , 1 = 1, NR ) 

READ (5,44) (ZMC( I ) ,1=1, NR) 

READ (5,44) (ZTC ( I ) , 1=1 , NR ) 

READ (5,44) (ZOC ( I ) , 1 = 1 ,NR ) 

WRITE (6,45) TITLE 

WRITE (6,46) ETA, LAMDA, 0P1,0P2, TNLMT 


DO 2 1 = 1, NR 

WRITE (6,47) I ,R I ( I ) ,R0( I),TI(I),TM(I) » T0( I) ,KIC( I ) ,KTC( I ) ,KOC( I ) , 
1Z MC( I ) , ZTC ( I ),ZOC( I > 

CALC. OF BLADE ELEMENT PARAMETERS 
E TA=ETA/DEGR AD 
L AMDA=L AMD A/DE GRAD 
DO 3 1 = 1, NR 
K IC(I) = KIC( I )/DEGRAD 
K TC ( 1 1 =KTC ( I )/DEGRAD 
KOC( I ) =KOC ( I )/DEGR AD 
E IC=0.0 
DO 4 1=1, NR 
S INKJC = SIN(KTC( I)) 

S INKOC=SIN(KOC( I ) ) 

S INK IC = SIN(K IC( I )) 

TNALP=( ROC I )-RI ( I) )/ZOC ( I ) 

TALP4 I )=TNALP 
ALP(.I)=ATAN(TNALP) 

CALPU )=S0RT(1./(TALP( I) **2+1. )) 

SALPtI )=TALP ( I ) *CALP 1 1 ) 

SNALJ?= SALP ( I ) 

C SAL P=CAL P ( I ) 

R ICm=RI( I ) /SNALP 

R OC ( I ) =RIC( I )+ZOC( I I/CSALP 

RTC( I)=RIC( I )+ZTC( I I/CSALP 

RMC( I)=RIC( I )+ZMC( I ) /CSALP 

C IC( X ) = ( S INK IC — SINKTC) /(RTCfl) — RIC <1)1 

COC( l)=(SINKTC-SINKOC)/(ROC(I)-RTC(I ) ) 

T 1=Z.TC( I) 

T2=ZMC( I) 
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on \j\ n * 


IF LTI.E0.T2) ICASE<I)=1 
IF (.TI.LT.T2) I CASE ( I ) =2 
IF IT1.GT.T2) ICASEC I )=3 

CALL RAEP IRIPI I),RISI I),EIPII),EISm .RICLI) .EIC,KIC< I),TIf I) ) 
BICI I)=CICC I )*RTC( I )+SINKTC 
BOCI J)=COC( I )*RTCI D+SINKTC 

ETC! I)=SUBF(KTC( I) ,KIC ( I) tBICL I),RTCU> .RICH) ) 

E0C(I)=ETCI I )+SUBF<KOC(I) tKTCII) tBOCII) »ROC( I ) tRTClin 
CALL RAEP (RDP( I ),ROS( I) ,EOP(I ItEOSII) »ROC 1 1 > »EOC( I ) *KOC( I it TO! I)) 
C 

CALL NEED (ICASE(I)tl) 

C 

CALL INTGR I It X INTI , XI NT2, X INT3 ) 

C 

T HECGI I )=X INT1/XINT3/SNALP 
CRCGI I )=X INT 2/XINT3 
RCGU)=CRCG( I)*SNALP 

T AE=TAN( ETA ) 

TNL=TAN(LAMDA) 

TAE2a=TAE**2 
T HECGO=THECG (NR ) 

ZCGOf(CRCGINR)-RICINR) )*CALP (NR) 

R CGO=?RCG( NR ) 

CALC. OF BLADE ELEMENT COORDINATES 
DO 6 I = it NR 
T 1=T AE 2 

T 2=RCG( I ) /RCGO 

T3=TAE/T2/( 1 .+T I )* ( SORT! T2**2* 1 1. + T1 )-Tl 1-1. » 

DELTil )=ARSIN(T3) 

Z CGI. I )=ZC GO + (RCGI I )-RCGO)*TNL 
T )=TALP( I ) 

T =S ALP C I) 

Z IC( I ) =ZCG( I )-( CRCGI I )— R IC ( I ) ) *CALP( I ) 

T4=DELT( I ) — THECGII ) 

THETA! I )=THECG0+T4 
XIC(.I)=RI( IJ*C0S(T4) 

X IPU)=RIP( I )*T*COS IE IP! I )/T+T4) 

X IS! I) =R I SI I )*T*C0SIEIS(I)/T+T4) 

XTC(I)=RTC( I >*T*COS(ETC(I )/T+T4) 

XTP(I)=RTPI I )*T*C0SIETPII)/T+T4) 

X TSI I ) = RTS( I )*T*COS( ETSI I ) /T+T4) 

XMCU)=RMC( I )*T*COS (EMC ( I )/T+T4) 

XMPI I )=RMP ( I )*T*COS(EMP( I ) /T+T 4) 

X MS I I )=RMS( I )*T*COS( EM SI I ) /T+T 4) 

XOC ( J )=ROC( I )*T*COS I EOC 1 1 ) /T+T 4) 

XOPI I )=ROP I I )*T*COS(EOP( I ) /T+T 4) 

XOSI I)=ROS( I )*T*COS( EOSI I ) /T+T 4) 

T 5=Z IC I Il-RICI l)*CALPII) 

T 6=CALPI I ) 

Z IP(J)=T5+RIP( I )*T6 
Z ISI I )=T5+R I SI I )*T6 
Z TCI.I ) = T5+RTC( I )*T 6 
ZTPf I) = T5+RTP( I ) *T 6 
ZTSl..I)=T5+RTS( I )*T6 
Z MCI I >=T5+RMC( I )*T6 
ZHPI J )=T5+RMP( I ) *T6 
Z MSI I )=T5+RM SI I )*T6 
Z OC 1 1 )=T5+ROCL I ) *T 6 
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Z OP ( I ) =T5+ROP( I)*T6 
Z OS LI )=T5+RO S( I >*T6 
C 

Y CG( I)=RCG( I )*SIN(DELT<n ) 

Y IC< I 1=RI ( I ) *SIN(T4) 

Y IPt I )=R I P( I )*T*SIN(EIP( I I/T+T4) 

Y IS! I)=RIS( I )*T*SIN(EIS(II/T+T4> 

YTC ( I)=RTC( I )* T*S IN (ETC ( I )/T+T4) 

Y TP ( U=RTP( I)*T*SIN(ETPm/TVT4l 
YT$m=RTS( I )*T*SIN<ETS(I)/T+T4> 

YHP ( 1 1=RMP ( I >*T*SIN(EMP(I)/T+T4) 

YMCm=RMC( I )*T^SIN(EHC( II/T+T4) 

YMSC1>=RMS( I )*T*SIN(EHS( I >/T+T4) 

Y OC( I ) =ROC ( I )*T*SIN(EOC(n/T+T4) 

Y OP ( I ) =ROP ( I )* T*SIN(EOP(I I/T+T4) 

YOSC I> = ROS( I )♦ T*SIN(EOS( I ) /T+T4) 

6 X CGI I)=RCG( I )*COS<DELT<I> ) 

C 

NX=NR 

no 7 K= 1* NX 
L^NRtl-K 

7 X(K)s?XCG(L) 

C 

DO 8 K= 1* N X 

CALL SINTP ( XIP *ZIP,NR*X(K) ,ZI PX(KM 
CALL SINTP < XIS*ZIS*NR*X(K) *ZISX(K) > 

CALL SINTP < XOP * ZOP * NR ? X ( K ) * ZO PX ( K ) ) 

CALL SINTP <XOS*ZOS*NR*X<K)*ZOSX<K)> 

8 CONTJNUE 
C 

CALL XM IN ( Z IP X * Z1 * NX ) 

CALL XH IN <ZISX*Z2*NX> 

Z MIN=AM IN l(ZItZ2) 

CALL XMAX (ZOSX*Zl*NX) 

CALL XHAX ( ZOPX*Z2* NX) 

Z MAX-AMAX 1 ( Z 1* Z2> 

DZ = ( ZM AX-ZM I N ) / 30« 

DZ= ADJ ( DZ ) 

Z X ( I 1=DZ*(AINT(ZMIN/DZ>— 1. ) 

DO 9 1 = 2* 32 

zx< n=zx( i-n+oz 

9 IF CZX( I) .GT.ZMAX) GO TO 10 

10 NZ= I 

c 

00 17 J= 1 * NZ 
DO 17 1=1 *NR 
T=SALPin 
T2=TALP( I ) 

RZ=RI( I l+( ZXCJ l-ZIC C 1 1 )*T2 
CAPRZ=RIC( I J4(ZX(J)-ZIC<I n/CALPII) 

IF (CAPRZ.GT.RTCdU GO TO 11 
SNCP=BIC( I )-CIC( I>*CAPRZ 
CALL FIX ( SNCP * C APPA ) 

EPC=SIJBF(CAPPA *KIC< I ) *B IC { I > ,CAPRZ,RICCI )) 

GO TO 12 

11 S NCP?80C{ I )-COC( I)*CAPRZ 
CALL FIX C SNCP * CAPPA 1 

FPC=ETC( I H-SUBFf CAPPA *KTC( I I ,BOC( I ) *CAPRZ,RTC(I ) ) 

12 IF ( CAPRZ •GT*RTP(I)I GO TO 13 
SNCP^BI Pi I )— CIPCI )*CAPRZ 
CALL FIX ( SNCP * CAPPA ) 
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EPP=E I P( I I + SUB F (CAPPA , KIP ( I ) »B IP (I ) ,C APRZ , RI P( I ) J 
GO TO 14 

13 SNCPttBOPI I )— COP( I )*CAPRZ 
CALL FIX ( SNCPtCAPPA) 

EPP=ETP( I ) + SUBF (CAPPA *KTP 11} »B(JP( I ) ,CAPRZ,RTP( 1)1 

14 IF (CAPRZ.GT.RTS(I) ) GO TO 15 
S NCP;=B I S( I)-CIS( I )*CAPRZ 
CALL FIX (SNCP, CAPPA) 

EPS=EI SCI )+SUBF( CAPPA, KIS( I ) *B I S( I ) »CAPRZ,RIS( I) ) 

GO TO 16 

15 SNCP=BOS( I)-CAS( D+CAPRZ 
CALL FIX (SNCP, CAPPA) 

EPS=ETS( I )+SUBF(CAPPA»KTS( I) ,BOS(I ) ,CAPRZ,RTS( I ) ) 

16 THETAC1 J, I )=THETA( I l+EPC/T 
T hFTAP ( J« I ) = THETA( I J+EPP/T 
T HETAS( J, I > = THETA( 1 )+FPS/T 

XP( J+I )=RZ*CCS(THETAP(J,I l-THECGO) 

XS( J,I )=RZ*COS(THETAS( J.I )-THECGO> 

YP(J,I >=RZ*SIN( THETAPI J,I l-THECGO) 

YS( J,1 )=RZ*SIN( THETAS! J, I l-THECGO) 

17 CONTINUE 
C 

C CALC. OF BLADE SECTION COORDINATES THRU BLADE ELEMENT 

C STACKING POINTS 

DO 20 K=1 »NX 
DO 19 J=1 »NZ 
DO 18 1=1, NR 

vm=xp(j,i) 

V 1 ( I l = YP ( J » I ) 

V 2 ( II = XS( J,I ) 

18 V 3( I J = YS ( J » I ) 

CALL SINTP ( V, Vl,NR,X(K) ,YPX( J) ) 

19 CALL SINTP ( V2,V3,NR,X(K) »YSX( J) ) 

WRITE (2) (YSX( J),J=1,NZ) »(YPX(J) »J=1,NZ) 

CALL SINTP ( ZX,YSX,NZ,ZISX(K) , YISX(K) ) 

CALL SINTP ( ZX, YSX , NZ, ZOSX( K) , YOSX(K) ) 

CALL SINTP ( ZX,YPX,NZ,ZIPX(K) , YIPX(K) ) 

CALL SINTP (ZX,YPX,NZ,ZOPX(K), YOPX(K)) 

20 CONTINUE 
REWIND 2 
TN0RM1=0. 

C 

C CALC. OF BLADE SECTION CENTER OF AREA 

DO 21 K= 1 , NX 

READ (2) I YSX( J ) » J=1»NZ) ,(YPXIJ)»J=1»NZ) 

21 CALL CGS (YCGX(K),ZCGX(K),YIPX(K),ZIPX(K),YISX(K),ZISX(K),YPX,YSX, 
1NZ, ZX, YOPX ( K ), ZQPX(K)»YOSX(K)» ZOSX ( K ) ) 

C 

C CALC. OF DIFFERENCES BETWEEN BLADE ELEMENT STACKING 

C POINTS AND BLADE SECTION CENTERS OF AREA 

DO 22 1=1, NR 

CALL SINTP ( X, ZCGX,NX»XCG( I ) ,ZCGR) 

CALL SINTP ( X, YCGX,NX, XCG ( I ) »YCGR) 

DRCG=( ZCGR-ZCGt I ) l*TALP( I ) 

XCGR=SQRT( ( RCG( I l+DRCG )**2-YCGR**2 ) 

CALL SINTP (X, ZCGX,NX, XCGR, ZCGR) 

CALL, SINTP (X,YCGX»NX, XCGR, YCG R) 

DRCG=( ZCGR-ZCG ( I ) >*TALP< I ) 
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X CGR»=SORT ( CRCGC I)+DRCG )**2-YCGR**2 ) 

V ( I )Sr7. CGR 

V 1( IA=YCGR 

V 21 I)=XCGR 

2? TNORMl=TNORMl+ABS(ZCGR-ZCG( I > ) +ABS ( YCGR— YCG (II ) 

WRITE (6,48) T NORM 1 , ( THECGt I ) * 1 = 1 » NR) 

WRITE (6,49) (CRCGII )»I=1»NR) 

IF (TNORMl.LT.TNLMT) GO TO 24 
C RE ALINEMENT OF BLADE ELEMENTS 

DO 23 I=L»NR 

THFCGI I)=THECG( I)-DELT(I )+ATAN(Vl( I)/V2(I) ) 

CRCGI I )=RIC( I1 + (VI I )-ZIC(I ))/CALP( I) 

23 RCGLJ )=SQRT( Vl( I )**2+V2( I )**2) 

REWIND 2 

GO TO 5 

24 REWIND 2 

IF (QP2.E0.0.) GO TO 27 

C PRINT OUT BLADE ELEMENT PARAMETERS AND COORDINATES 

DO 25 1=1, NR 
K IC( I )=KIC( I )*DEGRAD 
K IPI I)=KIP( I )*DEGRAD 
K ISC I )=KIS ( I )*DEGRAD 
KM( I »=KM( I )*CEGRAD 
KTC( I )=KTC{ 1 )*DEGR AD 
KTP( I )=KTP( I )*DEGRAD 
K TS ( I ) =KT S I I )*DEGRAD 
K 0C( II=KOC( I )*DEGRAD 
K OP ( I ) =KOP ( I )*DEGR AD 
K 0S( I ) =K0 S ( I )*DEGR AD 

25 ALP(I)=ALP( I )* DEGRAD 

WRITE ( 6, 50) (I , ALP ( I ) , KM ( I) , K IC ( 1 » »KTC ( I ) , KOC ( I ) , KI P( I ) ,KTP( I ) » KO 
IP ( I>»KIS( I ) ,KTS( I ) ,KOS( I ) ,I = 1,NR> 

WRITE (6,51 ) ( I ,C IC ( I ) ,COC ( I ) ,CIP( I ) ,COP( I) ,CIS( I) ,CAS( I ) , 1 = 1, NR) 
WRITE (6,52) 

DO 26 J = l» NZ 

26 WRITE (6,53) Z X ( J ) , ( I , YS ( J , I ) , XS ( J , I ) , YP ( J , I ) , XP ( J , I ) , 1= 1 , NR ) 

WRITE (6,54) 

write (6,55) (i,Yism,xis<n,zis(i),Yipm,xipm,zip(n,Yicm,x 
1IC( U,ZIC( 1 ) , I = 1 , NR > 

WRITE (6,56) 

WRITF (6,55) ( I,YMS( I) ,XMS( I) , ZMS(l) , YMP ( I ) , XMP ( I ) , ZMP ( I) , Y MC ( I ) , X 
1M C ( I ) , Z MC ( I ) , I = 1 »NR ) 

WRITE (6,57) 

WRITE (6,55) ( I ,YTS( I) ,XTS( I) ,ZTS( I) ,YTP(I) ,XTP( I) ,ZTP( I ),YTC( I ),X 
IT C ( I J , ZTC ( I ) » I = 1»NR ) 

WRITE (6,58) 

WRITE (6, 55) ( I , YOS ( I ) ,XOS( I) ,ZOS(I) , YOP ( I ) , X0P1 I ) ,ZOP( I) ,YOC( I),X 
10C1 I),ZOC( I ) , I = 1 »NR ) 

27 CONTINUE 

IF (0P1.E0.0.) GO TO 28 

C READ IN X-VALUES FOR BLADE SECTIONS 

NXQ=QP 1 

READ (5,44) <XQ(K) *K=1«NXQ) 

GO TO 30 

C BLADE SECTION X-VALUES AT BLADE ELEMENT STACKING POINTS 

28 NXO=NX 

DO 2.9 K= 1 , NX C 

29 X 0 ( K.1 = X ( K ) 

GO TO 34 

C CALC. OF UNROTATED BLADE SECTION COORDINATES 

30 DO 13 K=1,NX0 
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DO 32 J=1»NZ 
DO 31 1=1 , NR 
VlllffXPlJ, I > 

vi(u*vp(j,n 
V2i n=xs(Jti > 
v3( i*=ys« j,n 

CALL SINTP (V, Vl,NR,XO(K) ,YPX( JJ) 

32 CALL SINTP ( V2 , V3,NR ,XQ(K ) , YSX ( J) ) 

CALL. SINTP (XIP,ZIP,NR,XQ(K),ZIPX(KI ) 

CALL SINTP ( XIS,ZIS,NR,XQ(K> ,ZISX(K) ) 

CALL SINTP ( XDP ,ZOP,NR,XQ(K) * ZOPX(K» ) 

CALL SINTP ( XOS »ZOS,NR,XQ(K) , ZOSXI K) ) 

CALL SINTP (ZX,YSX,NZ,ZISX(K), YISX(K) ) 

CALL SINTP (ZX»YSX»NZ»ZOSX(K)» YOSX(K) ) 

CALL SINTP (ZX,YPX,NZ,ZIPX(K),YIPX(K) ) 

CALL SINTP ( ZX» YPX» NZ»ZOPX(K)» YOPX ( K) ) 

WRITE (2) (YSX(J),J=1,NZ),(YPX(J),J=1,NZ) 

33 CONTINUE 
REWIND 2 

C 

34 DO 35 K=1 • NXC 
ZICX(K)=.5*(ZISX(K)+ZIPX(K)) 

CALL SINTP (XMC,ZMC,NR,XQ(K),ZMCX(K)) 

CALL SINTP ( XMP »ZMP*NR»XO( K) »ZMPX( K) I 
CALL SINTP ( XMS,ZMS,NR,XQ(K),ZMSX(K) > 

ZOCXIK )=.5*( 2DSX(K)+Z0PX(K) ) 

CALL. SINTP (XTC»ZTC* NR »XO (K) , ZTCX ( K) ) 

CALL SINTP ( XTP »ZTP »NR»XQ(K)«ZTPX(K)) 

CALL SINTP ( XTS,ZTS,NR,XQ(K) ,ZTSX(K) ) 

Y ICXI K )=• 5* ( YI SX( K ) +YI PX ( K) » 

CALL SINTP ( XMC ,YHC ♦ NR,XQ( K) , YMCX( K> > 

READ (21 (YSX(J),J=1,NZ) ,(YPX( J).J=1,NZ> 

CALL SINTP ( ZX» YPXfNZ* ZMPX(K) » YMPX(K) ) 

CALL SINTP ( ZX,YSX,NZ,ZMSX(K), YMSX(K) ) 

YOCX(K>=.5*( YOSX(K)+YOPX(K) I 

CALL SINTP ( XTC,YTC,NR,XQ(K) ,YTCX(K) ) 

CALL SINTP ( ZX , YPX» NZ *ZTPX(K)» YTP X ( K) ) 

CALL SINTP ( ZX,YSX,NZ,ZTSX(K),YTSX(K)) 

JIX(K) = SORT( (ZISXIK)-ZIPX(K) )**2+ ( YI SX( KI-YI PXIK) 1**2) 

TOX(K)=SORT( (ZOSX(K)-ZOPX(K) )**2+ ( YOSX( K )— YOPX(K) ) **2) 

CALL SINTP ( XHC * TM . NR * XQ ( K ) ,TMX(K) ) 

CALL SINTP (XCG,YCG,NR,XQ(K),YCGX(K)) 

CALL SINTP ( XCG »ZCG , NR»XQ(K)» ZCGX ( K ) ) 

35 CONTINUE 
REWIND 2 

C PRINT OUT UNROTATED BLADE SECTION COORDINATES 

DO 36 K=l,NXO 

WRITE (6,59) XQ(K),ZICX(K), ZMC X(K) , ZTCX(K) »ZOCX( K) » ZI PX (K) • ZMPX ( K) 
1, ZTPXI K )» ZOPXIK ) » ZISX(K) » ZMSX( K) t ZTSX(K) » ZOSX ( K) , ZCGXI K ),YICX(K) ,Y 
2MCXI 1C) * YTCX ( K) » YOCX( K),YIPX(K) ,YMPX(K) , YTPX(K) , YOPX{ K) , YISX (K ) , YMS 
3X(K1,YTSX(K ),YOSX(K) ,YCGX(KI,TIX(K) ,TMX(K) ,TOX(K) 

READ (2) (YSX( J ),J=1,NZ) , ( YPX( J) , J=1 , NZ > 

DO 36 J=1,NZ 

36 WRITE (6,60) ZX( J) ,YPX( J I , YSX( J) 

REWIND 2 

C CALC. OF ROTATED BLADE SECTION COORDINATES 

DO 42 K=I,NXO 

READ (2) (YSX(J),J=1,NZ),(YPX(J),J=1,NZ) 

T l=ZOCX (K)-ZICX(K) 

T 2=.5*(TOX(K )-TIX(K) ) 

T3=Y0CX(K )— YICX(K) 
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SGX=I— T1*T2+T3*SQRT ( T 1**2+ T 3**2— T2**2 1 1 / (Tl**2+T3**2) 
G X=ARSIN( SGX » 

CGX-COS(GX) 

GAMXf K )=DEGR AD*GX 
T2=TIXIK>/2. 

DO 32 J=1,NZ 
T=ZXIJ)-ZICXIK» 

Vl( J» = IYSXt J l-YICXIK) >*CGX-T*SGX+T2 
V2(JI=IYPXIJ I— YICXIK > I *CGX— T* SGX+T2 
V 3( J )= ( YSX( J )— YICXIKJ J*SGX+T*CGX+T2 

37 V4{ JI=(YPX( J )— YICX(K) ) * SGX+ T+CGX+T2 
X HICfT 2 

XH0C*T0X(KI/2. 

T3=Y.ICX<K) 

T4=ZICXCK ) 

XHMC=*t YMCX(K)-T3)*CGX-CZMCXIK»-T4)*SGX+T2 
XHTC-I YTCXIK >- T3)*CGX- < ZTCXI Kl -T4) *SGX+T2 
XHIP*CYIPX(K )— T3 )*CGX— (ZIPX(K) -T4)*SGX+T2 
XHMP=< YMPX(K)-T3)*CGX-I2MPX(KI-T4)*SGX+T2 
X HOP*( YOPX(K)— T3I *CGX— ( ZOPX ( K )— T4) *SGX+T2 
XHTP*( YTPXIK )-T3)*CGX-IZTPX<K)-T4>*SGX+T2 
XHISsU YISXIK)-T3)*CGX-(ZI SX ( K >-T4 J +SGX+T2 
X HMS.= ( YMSXI K )-T3)*CGX- (ZMSX(K) — T4) *SGX+T2 
X HOS?=( YOSXIK )-T3)*CGX-(ZOSX(K)-T4)*SGX+T2 
X HT Sj=( YTSX(K >- T3 >*CGX- < ZTSX < K> -T4 ) * SGX+T2 
XHSP=t YCGCNR )- T3 >*CGX- ( ZCG < NR ) -T4) *SGX+T2 
XHCG=< YCGX ( K )-T 3 1 *CGX- < ZCGX < K > -T4 > *SGX+T2 
X L I G=T 2 

X LMCs=( YMCX ( K )-T3)*SGXUZMCX<K)-T4) *CGX+T2 
XLOC*=< YOC X( K )- T3 >*SGX+ ( ZOCX ( K) -T4) *CGX+T2 
XLTO( YTCX ( K )— T3)*SGX+IZTCX( K) — T4) *CGX+ T2 
XLIP=C YIPX(K>-T3)*SGX+(ZIPX(K>-T4) *CGX+T2 
XLMP=< YHP X ( K )— T3 )*SG X+ ( ZMPX ( K) -T4I *CGX+T2 
X LOP* I YOP X ( K J— T 3 )*SGX+ ( ZOPX C K) — T4) *CGX+T2 
XLTP-I YTPX ( K )— T 3 )*SGX + ( ZTPXIIO-T4) *CGX+T2 
XLIS»(YISX(K )— T 3 )*SGX+ ( Z1 SX I K ) — T4) *CGX+T2 
XLMS=( YMSXIK )-T3 )*SGX+ ( ZMSX ( K ) -T4) *CGX+ T2 
XLOS;=( YOSXIK )-T3)*SGX+(ZOSX(K)-T4)*CGX+T2 
XLTS*U YTSXI K )- T3)*SGX+ I ZTSX ( K ) -T4 ) *CGX+T2 
XLSPU YCGINR >-T3)*SGX+(ZCG(NR)-T4) *CGX+T2 
XLCG.-1 YCGX C K )-T3)*SGX + ( ZCGX < K ) -T4 > *CGX+T2 
C 

X L( 11=0. 

DL=X1.0C/53 . 

DL = ADJ(OL J 
C 

xpscn=XHic 

XHPI1>=XHSI 1 > 

X L ( 2 J = XL I C 
XU 3I = DL 
DO 38 J=4» 55 
L OUT*J 

X L l J l=XL( J-l M-OL 

38 IF (XL< J1 .GE.XLOC) GO TO 39 

39 L OUT—LOUT + 1 

XU LOUT- 1 )=XLOC 
X L( LOUT ) = XLOC+XHOC 
X HS (jLOUT ) = XHOC 
X HPI.LOUT) =XHS( LOUT I 
JM=LOUT-l 
DO 40 J=2» JM 
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CALL SINTP ( V3» V1,NZ,XL(J) »XHS< J) > 

40 CALL SINTP ( V4, V2, NZ,XL ( J » ,XHP ( J) ) 

C 

CALL CGS1 
BETA«=BETA*DEGRAD 

C PRINT OUT ROTATED BLADE SECTION COORDINATES 

WRITE (6,61) X0( K ) , GAMX( K ),TIX(K) , XLSP , LB AR, ARE A, I HI N, ILLCG , PHLCG, 
1ILL,J?HL,TMX(K) ,TOX( K ) , XHSP.HBAR, BETA, IMAX, IHHCG, IHH,XL IC, XLMC, XLTC 
2, XL DC, XLIP,XLMP,XLTP» XLOP, XLI S* XLMS, XLTS , XLOS , XLCG,XH IC,XHMC, XHTC, 
3X HOC, XH IP , X HMP , XHTP , XHOP , XH I S, XHMS , XHTS , XH OS , XHCG 
DO 41 J = 1 , LOUT 

41 WRITE (6,60) XL ( J ) , XHP ( J ) , XHS( J) 

42 CONTINUE 

C 

GO TO 1 

C 

C 

43 FORMAT (12A6) 

44 FORMAT (8F10.5) 

45 FORMAT ( IH1//12A6) 

46 FORMAT ( 1HK/30X, 34H INPUT FOR BLADE COORDINATE PR0GRAM//34X,3HETA,7 

1X.5HLAMDA, 5X ,3H0P1 , 7X, 3HOP2 , 7X , 5HTNLMT /29X, 5F1 0. 5/ /I X,7HEL EMENT, 4X 
2, 2HRI, 9X, 2HR0, 9X,2HTI ,9X» 2HTM, 9X,2HT0,9X,3HKIC ,8X, 3HKTC * 8X » 3HK0C ,8 
3X , 3HZMC* 8X,3HZTC,8X, 3HZ0C ) 

47 FORMAT ( 3X , 1 2, 2X , 11F 11 . 5 ) 

48 FORMAT ( IHJ/43H BLADE ELEMENT STACKING PARAMETER — TN0RM1 =,G10.3// 
16 H T-HECG/8G16 . 7/8G16.7/8G16.7) 

49 FORMAT (1HJ/5H CRCG/ 8G16. 7/ 8GI 6. 7/ 8G16. 7 J 

50 FORMAT ( 1H1/ /20X.20HBLADE ELEMENT ANGLE S//1X, 7HELEMENT ,4X, 3HAL P, 8X 
1, 2HKM* 9X» 3HK IC , 8X, 3HKTC , 8 X , 3HK0C , 8X ,3HKI P,8X , 3HKTP,8X,3HK0P, 8X, 3HK 
21 S, 8X» 3HKTS, 8X,3HKOS/( 3X,I2,2X,11F11.5)) 

51 FORMAT { 1HK , 20X , 24HBLADE ELEMENT CURVATURES//1X,7HELEMENT ,4X,3HCIC 
l, 8X, 3HC0C , 8X * 3HCI P, 8X,3HC0P, 8X ,3HC I S, 8 X , 3HC AS/ (3X, I2,2X,6F11.5) ) 

52 FORMAT ( 1H1/ /20X , 25HBL ADE ELEMENT COORDINATES) 

53 FORMAT (1HL/30X,3HZ =, F 10 . 5 , // 20X, 7HELEMENT, 1 5 X, 2H YS , 12 X, 2 HXS , 28X, 
12HYP,12X,2HXP/( 22X , I 2 , 10X , 2F14 .5 ,16X , 2F 1 4. 5 ) ) 

54 FORMAT ( 1HL / IX ♦ 7HEL EMENT , 9X, 3H YI S , 9X ,3HX I S ,9X , 3HZ I S , 13X, 3HY I P, 9X , 3 
1HXIP.,9X,3HZIP, 1 2X, 3HYIC , 9X, 3HXIC *9X,3HZIC) 

55 FORMAT ( 3X, I 2, 6X , 3F 12. 4, 4X , 3F1 2.4 , 4X , 3F1 2 .4 ) 

56 FORMAT ( 1HL / IX , 7HEL EMENT, 9X ,3H YMS, 9X,3HX MS ,9X, 3HZMS , 13X, 3HY MP, 9X , 3 
1HXMP*9X, 3HZMP, 13X, 3HYMC, 9X , 3HXMC , 9X ,3HZMC ) 

57 FORMAT ( 1 HL/ IX , 7HEL EMENT, 9X , 3H YTS , 9X, 3HX TS ,9X , 3HZTS , 13 X, 3HYT P, 9X ,3 
l.HXTP,9X, 3HZTP, 13X, 3HYTC , 9X, 3HXTC ,9X , 3HZTC ) 

58 FORMAT ( 1HL / IX , 7HEL EMENT, 9X, 3H YOS , 9X,3HX0S ,9 X , 3HZ0S , 13X, 3HY0P, 9X , 3 
1HX0P,9X,3HZ0P, 13X, 3HY0C,9X, 3HX0C, 9X , 3HZ0C) 

59 FORMAT ( 1H1/40X ,44HBLADE SECTION COORDINATES (UNROTATED) AT X = , F9 
1. 4/5X, 3HZ I C, 7X , 3HZMC,7X, 3HZTC , 7X, 3HZ0C ,7X,3HZIP,7X,3HZMP,7X,3HZTP, 
27X.3HZ0P, 7X, 3HZI S, 7X , 3HZMS, 7X, 3HZTS,7X , 3HZ0S ,7X,3HZCG/13F10 .4/5X , 3 
3HYIC, 7X* 3HYMC, 7X , 3HYTC , 7X , 3HY0C , 7X ,3HYI P .7X, 3HYMP, 7X . 3HYTP, 7X , 3HY0 
4P,7X.,3HYIS, 7X, 3HYM S, 7X, 3H YTS ,7X ,3H YOS , 7X , 3HYCG/13F 10 .4/5X, 2HT I , 8X, 
52HTM,8X, 2HTO/3F10.4/44X, 1HZ,9X,2HYP,8X,2HYS) 

60 FORMAT (39X.3F10.4) 

61 FORMAT ( 1H1/40X,42HBLADE SECTION COORDINATES (ROTATED) AT X =,F9.4 
1/ 5X,.5HGAMMA , 5X , 2HT 1 , 8X,5HL ( SP) ,5X, 5HL-BAR, 5X,4HAREA,8X ,4HI MIN, 8X, 5 
2H ILLCG, 7X, 5HPHLCG, 7X ,5HI (LL ) ,7X,3H PHL/4F 10.4,2 X,6G 12.4/ 5X ,2HTM,8X, 
32HTO, 8X,5HH( SP ) , 5X, 5HH-BAR , 5X, 4HB ETA, 8X , 4HI MAX , 8 X, 5H IHHCG, 19X , 5HI( 
4HH1/AF 10. 4, 2X» 3G12 .4, 12X, G1 2.4/5X, 5HL ( I C ) , 5X , 5HL ( MC ) , 5X , 5HL ( T C ) , 5X 
5, 5HL1 OC ) * 5X » 5HL ( IP ) , 5X , 5HL ( MP) ,5X, 5HL ( TP ) , 5X , 5HL ( OP ) , ?X , 5HL ( I S ) , 5X 
6, 5HL4 MS ), 5X , 5HL ( TS ) , 5X, 5HL ( OS) ,5X, 5HL ( CG ) /13F 10.4/ 5X ,5HH( IC),5X,5H 
7H ( MC I , 5X, 5HH (TC ) ,5X, 5HH( OC ) , 5X , 5HH (IP),5X,5HH(MP),5X,5HH(TP),5X,5H 
8H ( OPJ , 5X« 5HH( IS) ,5X,5HH(MS) ,5X,5HH(TS) ,5X,5HH( OS) ,5X,5HH(CG)/13F10 
9. 4/4AX , 1HL , 9X, 2HHP, 8X, 2HHS) 

END 
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SUBROUTINE ITER < K2 ,C , B , K l , E 1 * Rl ♦ E2 f R2 , XK ) 

REAL Kl ? K2 

Z=E2-E1 

K 2* XX 

IGO-J. 

I T=0 
DK= • 1 

1 C = (SINtKl) — SINIK2I )/ (R2— R 1 ) 

B-C*Rl+SI Nf K 1 ) 

F*Z— &UBFI K2fKl»B»R2»Rl) 

I T«IJ>1 

IF f IT «GE .101 IG0=3 
GO TO I 2f 3t 6 ) « IGO 

2 T K=ICi2 
TF~F 
K2-TJC+DK 
I GO- 2 
GO TO 1 

3 IF <F*TF.LE.O.O) GO TO 5 

IF CABSm.GT.ABS(TF)) GO TO 4 
TK=K^ 

TF*F 

K2=TK+DK 
GO TO l 

4 OK=-OK 

K 2=TK + DK 
GO TO 1 

5 IF <ABS(R2*F).LT..0001) IG0=3 

DK~TF*OK/(TF-F ) 

K ?=TK+DK 
GO TO 1 

6 RETURN 
END 


SUBROUTINE I TER I (KT*RTfET» KH» RM* E N f B »C» KTC t RTC * ETC! 
REAL, KM»KT « KTC 
I G0*1 
I T=0 
RT^RTC 

1 SNK TASINI KM)+C*(RM-RT> 

CALL FIX (SNKT*KT) 

ET=EN+SUBF<XT t KH 9 Bt RTtRM) 

R TC2jpRT /EXP ( TAN I KTC I * ( ETC-ET ) I 
F =RT£— RTC2 
I T= I 1+ 1 

IF UT.GE.10) I G0=3 
GO TO t 2v 3* 6 )* IGO 

2 R2=RJ 
F2=F 
D EL*F 
RT*RJ+OEL 
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1 60 =£ 

60 Tit 1 

IF IR*F2.LE.0.0> GO TO 5 
IF <ABS(FI.GT.ABS(F2>> GO TO 4 
R2=RT 
F2=F 

RT=RT+DEL 
GO TO 1 

4 0£L*^DEL 
R T=R£+DEL 
GO TO 1 

5 IF U8S( F2I.LT. .0001 1 IG0=3 

D£L»*2*DEL/{ F2-F) 

RT=R«+DEL 

GO TO 1 

6 RETURN 
END 


SUBROUTINE SINTP C Z, W,N,X1,Y1I 
DIMENSION X ( 56 ) • YI56I, ZI56), W(56) 

00 1 1=1, N 

xm*z<n 

1 yinxwm 

CALL 50RTXY IX,Y,N) 

C 

K *1 

DO 2 I=l,N 

IF (il.GT.Xdll GO TO 2 
IF MU.EQ.XI II I GO TO 3 
JF LX1.LT.XC I)) GO TO 4 

2 K»I *4 
GO TO A 

3 Y1=Y4K) 

RETURN 

4 IF IJ^.EO.1) GO TO 5 

IF Ift.GT.N) GO TO 6 

IF ifc.EQ.N) K=N— 1 

W 1=(J<1-X(K ) J*(X1-X(K+W)/(XIK-1)-X(KJ ) /{ XTK-1 J -XIK+1TI 
N2=(J(l— X( K— 1 ) 1*1 XI— XIK+1) ) / ( XI XI— X IK— II I /{ XI K) — XI K+l ) ) 
N3*i^i-x(K-in*(xi-x(Kn/<xiK+i»-x(K-in/{xiK+n-xnm 

Yl=YiK-l>*Hl4YCKI*N2+Y(K+ll*M3 

RETURN 

5 J*1 
Lf*2 

GO TO T 

6 J=N-1 
L*N 

7 Y 1=YJ,J 1+1X1— X(JI1*IYCLI— YCJ)I/(XIL| — XIJI1 
RETURN 

END 
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SUBROUTINE CGS ( YCG *ZCG *YI PX*Z I PX» YI SX, ZI SX *YPX* YSX* NZ * ZX* YOP* . ZOP 
IX « YOSX 9 ZO SX I 

c 

0 I MENS I ON V C 51 ) » VU51), V2I51), V3(51), YPXI50) * YSXI50)* ZXC50) 

C 

DYI=OriSX-YIPX 
DZI=*ISX-ZIPX 
DZO=ZQSX— ZOP X 
OYO=vYOSX-YOPX 
C 

T I2=DY I**2+DZI **2 
T02=DY0**2+DZO**2 
C 

T IX*iSQRT( T I 2 I 
TOX=£ORT( TO 2 I 
C 

SINB<I = -DZI/TIX 
S INB0=-DZ07T0X 
COSBJ=OYI/TIX 
C OS BO= 0 YO / TO X 
C 

K Z 1^2 

IF tOZI.GE.O.O) KZ I = 1 
C 

GO TO ( 1«2)«KZI 

1 ZHIN*ZIPX 

Z M !N2=Z IS X 
GO TO 3 

2 Z*Ifi*ZISX 
ZMIN2*ZIPX 

3 K ZO=2 

IF (OZO.GE.O.O) KZO=l 
C 

GO TO ( 4, 5 ) * KZO 

4 ZMAX*ZOSX 
ZMAX,2=Z0PX 
GO TO 6 

5 ZMAX^ZOPX 
ZKAX2=ZOSX 

C 

6 DZ=UMAX-ZMIN) /50. 

Z=ZHIN 

C 

DO L4 1=1*51 
IF U.GE.ZMIN2) GO TO 9 
GO TO ( 7* 8 1 * KZ I 

7 Y S=YJPX+C Z-Z IPX )*DYI /DZI 
CALL SINTP ( ZX, YPX»NZ*Z* YP) 

GO TO 13 

8 Y P=YrISX+C Z— Z 1SX )*DYI /OZI 
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CALL SINTP < ZX, YSX,NZ*Z,YS) 
GO TO 13 

IF ( Z *GT « ZMAX 2 ) GO TO 10 
CALL SINTP I ZX,YSX,NZ,Z*YS> 
CALL SINTP C ZX.YPX,NZ*Z,YP) 
GO TjOl 13 

GO TjQ ( 11 1 12 )f KZG 
YP=YGPX+CZ-ZOPXI*OYO/DZO 
CALL SINTP (ZX,YSX,NZ.Z.YS 1 
GO TO 13 

YS=Y£SX4( Z-ZOSX)*DYO/OZO 
CALL SINTP CZX^YPXtNZ,Zf YP) 

V(I^YS-YP 

Vl( U=,5*(YS4YP)*V< I > 

V 2( 1J=Z*V( I ) 
t= Z+DZ 

PIRA 0 = 3. 1415927 

T*PI$AD/ 8.0 
A I*T 42 *T 
A 0 =TQ 2 *T 


T*2~/3./PIRAD 

YBI»4YISX4YIPXI /2.-T*TIX*SIN8I 

Y B0^1Y0SX4Y0PX ) /2.4T*T0X*SINB0 

ZBI=4ZISX+ZIPX)/2.-TIX*COSBI*T 
Z B0=4Z0SX+Z0PX ) /2# + T*TOX*COSBO 

CALL FNTGRL <51.DZ,V,V3I 
A =V31 51 )4A 1 4 AO 

CALL FNTGRL (51 ,DZt VI * V3 > 

Y CG=1 V 3( 5 1 ) 4 YB I* A I4YBO+A0 1 /A 

CALL FNTGRL C51 * OZ t V2 f V3 ) 
ZCG=4V3(51 14ZB I*AI4ZB0*A0 > /A 

RETURN 

END 


SUBROUTINE FIX IA f B) 

IF CA8S(A).LT.l*) GO TO 1 

IF U.LT.O*) T~-l. 

B=T*1. 570795 

R ETURN 

8=ARSIN(AI 

RETURN 

END 



SUBROUTINE C6SI 

DIMENSION X I 56 ) » HSI56), HP(56I, VI56J, Vl<56), V2(56), V3(56), XI 
11561 

COMMON /EXTRA/ X»HS» HP »N»LBAR* HBAR , ILL* IHH ,PHL» AREA, ILLCG, IHHCG, PH 
ILCG.BETA, IM IN* IMAX 

REAL LBAR, ILL* IHH, ILLCG, I HHCG, I MIN, I MAX 
DX=XXN>/49. 

Xl( u=o. 

00 1 1 = 1,50 

CALL SINTP <X,HS,N,Xl(I»,Vm» 

CALL SINTP (X,HP, N, Xl(I), VIII) I 

1 XIU«I=X1II)+0X 
DO 2, 1 = 1, 50 

2 V2CU=VI II-VK I > 

C ALL, FNTGRL I50,DX,V2, V3I» 

ARFAsfV3I 50) 

DO 3. 1 = 1,50 

3 V2( U=V2I 1 1*X1 Cl) 

CALL FNTGRL 150, OX, V2, V3 I 
L BAR»V3( 50 ) /AREA 
DO 4 1=1,50 

4 V2( IJ=V2( I)*X1|I) 

CALL. FNTGRL ( 50 , DX , V2 » V3 » 

1 HH=.V3 ( 50 I 
DO 5 1 = 1, 50 

5 V2I I.4=V( I I**2-V1( I 1**2 
CALL FNTGRL I 50, DX , V2 , V3 > 

HBAR*V3( 501/2. /AREA 

DO 6 1=1, 50 

6 V 2( IJ = V2( I l*Xl C I ) 

CALL, FNTGRL 150 ,DX, V2, V3 > 

PHL=W3( 501/2. 

DO 7 1=1, 50 

7 V2i n=vm**3-vim**3 
CALL FNTGRL < 50 , DX , V2, V3 > 

I LL=V 3< 501/3. 

I HHCG= IHH-AR EA*LBAR**2 
I LLCG= ILL -ARE A* HBAR** 2 
PHLCG=PHL-AREA*HBAR*LBAR 
B FTAsrAT AN C 2.*PHLCG/( IHHCG— ILLCG1 ) 

IMIN*I ILL CG+ IHHCG ) /2.-M I LLCG-I HHCG ) /2. *COS< BET A ) -PHLCG*S I N< BET A ) 

I MAXa=C ILLCG* IHHCG)— I MIN 
BETA*BETA/2. 

R ETURN 
END 


SUBROUTINE XMAX (X,XM,N) 
DIMENSION X(IOO) 

XM=XU I 
DO L 1 = 2, N 

IF (Xm.LT.XM) GO TO 1 
XM=X«I) 

1 CONTINUE 

R ETURN 
END 
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SUBROUTINE XMIN IX.XM.N) 
DIMENSION X(IOO) 

XM=X.I 1 ) 

DO 1 I*2.N 

IF IXID.GT.XM) GO TO 1 
XM=)Uf ) 

1 CONTINUE 

RETURN 
END 


SUBROUTINE NEED IIC.I) 

DIMENSION ALPI24), B0CI24), BICI24), BIPI24), BISI24), B0PI24), BO 
IS ( 241* CISI24), COPI24), CASI24), CICI24), EIPI24), EISI24), EMCC2 
24), SMPI24), EMSI24) , EOCI24), E0PI24), EOSI24), ETC(24), ETPI24), 
3 ETSX24), KICI24), KIPC24), KISC24), KMC24) , KOCI24), K0PC24), KOS 
4I24U KTC ( 24 ), KTPI24), KTSI24), RI(24I, R0(24) , RIC(24), RIPI24), 
5 RISI24), RMCI24), RMPI24), RMSI24), ROCI24), R0P(24), R0SI24), RT 
6CI24I, RTPI24), RTSI24), TII24), TMI24), T0I24), ZMCI24), Z0CI24), 
7 ZTC124), C0C(24), CIPI24) 

COMMON AL P, BOC ,BIC»BIP,BI S , BOP , BOS ,CI P »CIS»COP»CAS»CIC,EIP»EIS«EMC 
1, EMP.EMS, EOC,EOP,EOS,ETC ,ETP ,E TS, KIC ,KI P ,K I S ,K M, KOC, KOP , KOS ,KTC, KT 
2P ,KT.6* RI,RO,RIC,RIP,RIS,RMC»RMP»RMS»ROC, ROP,ROS*RTC , RTP ,RTS , T I , TM» 
3T O, ZMC , ZOC * Z TC , NR, COC 

REAL KIC, KIP, K IS, KM, KOC , KOP ,K0 S, KTC ,KTP ,KTS , LA MO A 
C 

C CALC. OF BLADE ELEMENT SURFACE PARAMETERS 

C 

GO TO 11,2,3), IC 
C 

C TRANSITION AT MAX, THICKNESS 

C 

1 KMIU=KTCCI) 

EMC ( J ) = ETC I I I 

CALL RAEP ( R MP ( I ) ,RMSf 1 1 , EMP ( I ) ,EMSI I ) ,RMC II) , EMCCI) ,KMII) ,TM! I)) 
KTPf J )=KTC( I ) 

K TSC I )=KTC( I ) 

ETPU) = EMPI I ) 

ETSU) = EMS( I ) 

RTPIJ)=RMP< I ) 

RTSI I ) *RM SI I ) 

C 

CALL ITER IKIPI I),CIP( I) ,BIPtI ),KTP(I ) ,ETPCI) ,RTP( I) *EIP( I ),RIP( I) 

1 « K I C4- 1 ) ) 

CALL ITER (KISI I),CISm,BIS(I I.KTSII) ,ETStI) ,RTS(I),EISm ,RIS( I) 
1, KICI I ) ) 

CALL ITER CKOPC I ) .COP I I ) ,BOP ( I ) ,KTP(I ) , ETP(I) ,RTP{ I) ,EOP!I ) ,ROP(I) 

I , KOC.l I ) ) 

CALL ITER ( KOS( I ) ,CAS( I ) ,BOS( I ) , KTS ( I ) ,ETS ( I ) , RTS(I) , EOS (I) ,ROS ( I ) 
1.K0C4 I) ) 

R ETURN 
C 

C TRANSITION AHEAD OF MAX. THICKNESS 

C 

2 S INKM=BOC( I )-COC( I )*RMCI I ) 

K Ml IJ =ARS INI SINKM) 

EMC C J ) = ETC I I )+SUBFIKMII),KTCII ) ,BOCII ) ,RMCII) ,RTCII) ) 

CALL RAEP CRMPI I > ,RMSI I) , EMP 1 1 ) ,EMS< I ) , RMC 1 1 ) , EMC 1 1 ) ,KM{ I ) , TMI I)) 
CALL ITER I KOP I I ) ,COP{ I ) ,BOP 1 1 >,KMtI),EMP(I ) ,RMP 1 1 ) ,EOP 1 1 ) , ROPII), 
1K0CIJ) ) 

CALL ITER IKOSI I) ,CASI I ) , BOS 1 1 ) , KM 1 1 ) ,EMS ( I ) ,RMS 1 1 ) , EOS 1 1 ) , ROS 1 1 ), 
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ui n in o 


IKOCC I) » 

CALL I TER 1 CKTPC II .RTPCII ,ETPC I I ,KMII 1 , RMPC I 1 r EMP( 1 1 ,BOPC I 1 ,COPC II 
l.KTCU l.RTCC II ,ETCC III 

CALL ITER 1 CKTSC Il»RTSCIi,ETSC I). KM II) .RHSCII ,EMS(II ,BOS(I) ,CAS( II 
1,KTC4I!,RTCC IltETCIIll 

CALL ITER (KIP( IltCIPCII.BIPCI ) »KTPC I I ,E TPCI 1 ,RTP C 1 1 ,E IPCI I , RI PC II 
I* Kid III 

CALL ITER (KISI II ,C I S( II ,B I S II I ,KTSC 1 1 ,ETSII I , RTSC II ,E ISCI) .RISC II 
L.KICCI 1 1 
R ETURN 

TRANSITION BEHIND MAX. THICKNESS 

SINKM*BICC I)-CIC(I>*RMCm 
K HI I |* ARS INI SINKMI 

EMCC I)=ETCC I l+SUBFCKMC II.KTCCI1 .BICCI 1 ,RMC(I> »RTC( 1 1 I 
CALL RAEP CRMPC II.RMSCII.EMPCI l.EMSCI I ,RMCCI 1 , EMC Cl) ,KMCI I ,TM(I) I 
CALL ITER (K!P(I)*CIP(I)tBIPd )»KM(I I »EMP( II ,R MP (1 1 ,EI PC 1 1 , RIP CI1, 
IK ICC II I 

CALL ITER ( K IS C IltCISC Il.BISCI ) » KM ( I ) ,EMS( 1 1 .RMSCI I ,EIS (1 1 ,RIS ( I I. 
IK ICC II I 

CALL ITER 1 CKTPC II .RTPC I I ,ETPC 1 1 .KMCI I .RMPC II .EMPCIl ,BIPC 1 1 ,CIPC II 
ltKTCCII.RTCC I) »ETCC 111 

CALL ITER 1 CKTSC II .RTSC 1 1 ,ETS< 1 1 , KMC 1 1 ,RMS C 1 1 , EMS C 1 1 ,B IS C 1 1 ,C IS CII 
l.KTCil l.RTCC 1 1 * ETC €11! 

CALL ITER CKOPC 1 1 ,COP C 1 1 . BOP Cl l.KTPCII .E TP Cl I ,RTPC 1 1 .EOPCI 1 ,ROPC II 

i , kocii 1 1 

CALL ITER CKOSC I I.CASC 1 1 .BOSCI 1 .KTSCII .ETSCI1 , RTS C 1 1 .EOS C 1 1 ,ROSC 11 
l * KOCi I I 1 
RETURN 
END 


FUNCTION SUBF (X.XO.B.R.ROl 
IF CX.EO.XOl GO TO 3 
B?=B*B 

IF CB2.LT. 1. 1 GO TO I 
IF (B2.GT.1. 1 GO TO 2 
T = 1 . 

IF LB. EO.-l . 0) T =- 1 . 

T l=.7854+T*X/2. 

T2=.T854+T*XO/2. 

SUBF*CX-X01-T*C TANCT11-TANCT2) 1 
RETURN 

1 T =T ANC X/2 • I 

T l=TANCXO/2. 1 
T 2= SORT ( 1 .— B 21 

SUBF*C X-XO I +B/T2*C ALOG ( ABSC C-B*T+1.-T21/C-B*T+1.+T2I II-ALOGCABSC (- 
1B*T 1+1 .-T21 / C-B*T1+1.+T2 1 111 
RETURN 

2 T =T ANC X/2 • 1 

T 1=TAN C XO/2. I 
T 2* SORT C B2- 1 .1 

SUBF=2.*B/T2*C ATANC l-B*T+ 1 . 1 /T 21-ATAN C C - B*T1+1 . 1 /T21 1+CX-XOl 
RETURN 

3 SUBF*TANCXI*ALOGCR/RO! 

RETURN 

END 
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SUBROUTINE INT6R (L,X1,X2,X3) 

DIMENSION ALP(24), B0C(24), BIC(24), BIP(241, BIS(24)» B0P(24), BO 
IS(24*. Cl S( 24) * C0PI24), CAS<24), CICI24), EIP124), EIS(24), EMC (2 
24), EMP(24), EMSI24), EOCI24), E0PI24), E0SI24), ETC124), ETP(24), 
3 ETS.I24), KIC( 24), KIPI24), KISI24), KMI24), K0CC24), KOPI24), KOS 
4124).* KTC ( 24 ), KTPC24) , KTSI24), RII24), R0(24) , RICI24), RIPI24), 
5 RIS124), RMC124), RMPI24) , RMSI24), R0CI24), R0PI24), R0SI24), RT 
6CI244, RTPI24), RTSI24), TII24), TM124) , T0I24), ZMC124), Z0CI24), 
T. ZTCA24) , COCI 24) , CIPI24) 

COMMON ALP, BOC, BIC, BIP ,BI S,BOP ,BOS »CI P»CI S»COP ,CAS,C IC,EIP,£IS,EMC 
1, EMP,EMS» EOC,EOP,EOS»ETC,ETP»ETS,KIC ,KI P »KIS»KM»KOC*KOP,KOS »KTC,KT 
2P ,KTS,RI.RO,RIC*RIP»RIS,RMC ,RMP,RMS,ROC , ROP, ROS,RTC, RTP.RTS, TI ,TM, 
3T0,ZNC,Z0C,ZTC,NR,C0C 

REAL KIC, KIP,KIS,KM,KOC,KOP , KOS, KTC ,KTP,KTS,LAMDA 
DIMENSION FlUOl), F2(10i), F3(10I), V(lOl) 

RMAX*ROC(L) 

RMINcRIC(L) 

DR=C.iRMAX-RMIN) /100. 

R-RMIN 

SNKTJB=SIN(KTS(L) ) 

SNKT.R=SIN(KTP(L ) ) 

C 

DO 5 1*1,101 

IF LR.GE.RTSCL) ) GO TO 1 

SNKS*BIS(L)-R*CIS(L) 

XKS*ARSIN( SNKS ) 

EMAX*ETS(L )+SUBF!XKS,KTS(L) ,BIS(L) ,R,RTS(L>) 

GO TO 2 
C 

1 SNKS*BOS( L )— R*CAS( L ) 

XKS*.AR SIN (SNKS) 

EMAX*ETS( L ) +SUBFI XKS,KTS CL) ,BOS(L) ,R,RTS(L) ) 

C 

2 IF (R.GE.RTP(L) ) GO TO 3 

SNK P*8IP( L )-R*C IP( L ) 

XKP*AR SIN ( SNKP ) 

EMINieETPI L ) + SUBF CXKP.KTP (L),BIP(L) ,R,RTPCL)) 

GO TO 4 
C 

3 SNKP*BOP(L )-R*COP(L) 

XKP*AR SIN (SNKP) 

EMIN*ETP( L )*SUBF(XKP»KTP(L) ,BOP(L) ,R,RTP(L) ) 

4 F 3( I )*C EMAX-EM IN )*R 
F2C I l=R*F3( I ) 

F1(II*,5*(EMAX+EMIN)*F3(I) 

^ 

CALL FNTGRL (101, DR, FI, V) 

X 1*V1101) 

CALL FNTGRL (101,DR,F2,V) 

X 2*VA101 ) 
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C ALL, t-FNTGRL 1101. DR, F3, VI 
X3=V<101> 

RETURN 

END 


FUNCTION ADJ (0) 

A*l«43 

1 IF Ut.GT. 1.01 GO TO 2 
A*A/10. 

D =D ^<1 0 • 0 
SO TO 1 

2 IF (O.LE. 10*01 GO TO 3 
A*A*40.0 

D=D/A0,0 
GO TO 2 

3 IF (A*0«LT*D -AND *D *LE* 2*0 ) ADJ=2.0*A 
IF (£*0*LT*D .AND.D *L E. 5*0 ) A0J=5.0*A 
IF tC.b.LT.D.AND.D.LE.lO.O) ADJ=10.0*A 
RETURN 

END 


SUBROUTINE RAEP (RP f RStEPtES »RC »EC tXKC* TC) 
IF iXKC + E Q*0 *0 ) GO TO I 
V =TC/2 .*S INI XKC ) 

R P=R£+V 
RS=R£-V 
V -COJAN( XKC ) 

EP=EC+V*ALOG(RC/RP J 
ES=EC+V*ALOG (RC/RS ) 

RETURN 
1 R P=RC 

R S=R€ 

V=T C/2 • /RC 

EP=E1>V 

ES=E£>V 

RETURN 

FNO 
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CT> 

00 

SAMPLE BLADE FOR BLADE COORDINATE PROGRAM 


INPUT FOR BLADE COORDINATE PROGRAM 

FTA LAMDA 0P1 0P2 TNLMT 

-1.00000 L. 00000 5.00000 1.00000 0.00010 


FLEMFNT 

R I 

RO 

TI 

TM 

TO 

1 

20.00000 

19.00000 

0.10000 

0.10000 

0. 10000 

2 

18. 00000 

17.80000 

0.10000 

0.10000 

0.10000 

3 

15.00000 

15.01000 

0.10000 

0.10000 

0. 10000 

4 

12.00000 

13.00000 

0.10000 

0.10000 

0. 10000 

5 

10.00000 

12.00000 

0.10000 

0.10000 

0. 10000 


BLADE ELEMENT STACKING P AR AMETER— TN0RM1 = 0.287E-01 
THECG 

0 C. 2764833E-02 0. 941 8329F-02 0 . 2063697E-0 1 

CRCG 

-43.605685 -179.B9311 3001.0398 27.955353 

BLADF FLFMENT STACKING PARAMETER — TN0RM1 = 0.-116F-02 
THECG 

0.297 1 57 IE- 04 0. 267641 IE-02 0. 9004 142F-02 0. 1988389F-01 

CRCG 

-43.604660 -179.89522 3001.0386 27.954099 

BLADE FLFMENT STACKING PARAMETER — TN0RM1 = 0.106E-03 
THECG 

0 .2260716E-04 0. 2676398E-02 0. 9004064F-02 0 . 1984415E-01 

CRCG 

-43.604774 -179.89527 3001.0386 27.954002 

BLADE ELFMENT STACKING PARAMETER — TN0RM1 = 0.407E-04 
THECG 

0.226 71 88F-04 C. 2676431E-02 0. 900401 8E-02 0. 1983933E-01 

CRCG 

-43.604771 -179.89527 3001.0386 27.953982 


BLADE ELEMENT ANGLES 

ELEMENT ALP KM KIC KTC KOC 

1 -26.56506 0. 0. 0. 0. 

2 -5.71060 0. 8.50000 0. -8.50000 

3 0.28648 0. 24.00000 0. -24.00000 

4 26.56506 0. 37.00000 0. -37.00000 

5 45.00002 0. 45.00000 0. -45.00000 


KTC KTC KOC ZMC ZTC ZOC 

0 . 0 . 0 . 1.00000 l . 00000 2.00000 

8.50000 0. -8.50000 1.00000 1.00000 2.00000 

24.00000 0. -24.00000 1.00000 l. 00000 2.00000 

37.00000 0. -37.00000 1.00000 1.00000 2.00000 

45.00000 0. -45.00000 1.00000 1.00000 2.00000 


0.3 865 755E-01 


15.570537 


0 • 38 24395F-01 


15.569048 


0 • 38232 25E-0 1 


15.568880 


0. 3823 1 1 7F-0 1 


15.568888 


KIP 

KTP 

KOP 

KIS 

KTS 

KOS 

-0.12886 

0 . 

-0.13253 

0.12886 

0 . 

0.13253 

8.46652 

0 . 

-8.53109 

8.53185 

0 . 

-8.46941 

23.99986 

0 . 

-24.00399 

23.99423 

0 . 

-24.00585 

37.20314 

0 . 

-36.80192 

36.80208 

0 . 

-37.19270 

45.36453 

0 . 

-44.65568 

44.64384 

0 . 

-45.33531 


Output Listing for Sample Blade 



> oJ in; 


FL FMENT 
1 


BLADE ELEMENT CURVATURES 


CIC 

COC 

CIP 

CUP 

CIS 

CAS 

0 . 

0 . 

-0.00201 

0.00207 

0.00201 

-0.00207 

0.14708 

0.14708 

0.14759 

0.14870 

0.14655 

0.14548 

0.40674 

0.40674 

0.41517 

0.41525 

0.39853 

0.39873 

0.53828 

0.53828 

0.55577 

0.55063 

0.52177 

0.52651 

0.50000 

0.50000 

0.51607 

0.50973 

0.48476 

0.49065 


BLADE ELEMENT COORDINATES 



l - -0.00000 




ELEMENT 

YS 

XS 

YP 

XP 

1 

-0.10329 

20.07893 

-0.20284 

20.07817 

2 

-0.13970 

18.01263 

-0.24102 

18.01156 

3 

-0.18675 

14.99844 

-0.29811 

14.99664 

4 

-0.22718 

11.98179 

-0.35505 

LI. 97868 

5 

-0.31272 

9.99511 

-0.45418 

9.98968 


ELEMENT 

7 = 

0.10000 

YS 

XS 

YP 

XP 

1 


-0.10276 

20.02893 

-0.20261 

20.02817 

2 


-0.12341 

18.00274 

-0.22456 

18.00177 

3 


-0.14173 

14.99943 

-0.25078 

14.99800 

4 


-0.14975 

12.03301 

-0.27015 

12.03091 

5 


-0.19367 

10.09814 

-0.32172 

10.09487 


ELEMENT 

Z = 

0.20000 

YS 

XS 

YP 

XP 

1 


-0.10225 

19.97893 

-0.20236 

19.97817 

2 


-0.10867 

17.99284 

-0.20965 

17.99195 

3 


-0.10180 

15.00025 

-0.20891 

14.99914 

4 


-0.08308 

12.08366 

-0.19779 

12.08232 

5 


-0.09354 

10.19957 

-0.21232 

10.19779 



l = 

0.30000 




ELEMENT 


YS 

XS 

YP 

XP 

1 


-0.10177 

19.92893 

-0.20208 

19.92817 

2 


-0.09545 

17.98291 

-0.19628 

17.98210 

3 


-0.06673 

15.00095 

-0.17217 

15.00011 

4 


-0.02594 

12.13391 

-0.13626 

12.13318 

5 


-0.00916 

10.30000 

-0.12132 

10.29929 



-3 

o 


7 = 0.40000 
YS 

0.10131 
0.08375 
0.03630 
0.02258 
0.06156 


2 = 0.50000 

YS 

0.10088 
0.07357 
C. 01033 
0.06314 
0.12006 


7 - 0.60000 
YS 

0.10047 
0.06489 
0.01128 
0.09625 
0.16732 


7 = 0.70000 
YS 

0.10009 
0.05770 
0.02869 
0.12228 
0.20405 


ELEMENT 

1 

2 

3 

4 

5 


FLEMENT 

1 

2 

3 

4 

5 


FLEMENT 

1 

2 

3 

4 

5 


ELEMENT 

1 

2 

3 

4 

5 


Z = 0.80000 


FLEMENT 

1 

2 

3 

4 

5 


YS 

-0.09972 

-0.05202 

0.04199 

0.14147 

0.23069 


XS 

19.87893 
17.97297 
15.00155 
12. 18392 
10.39982 


YP XP 

-0.20177 19.87817 

-0.18444 17.97222 

-0.14035 15.00094 

-0.08436 12.18365 

-0.04582 10.39990 


XS 

19.82893 

17.96302 

15.00209 

12.23378 

10.49931 


YP XP 

-0.20144 19.82817 

-0.17413 17.96232 

-0.11323 15.00167 

-0.04121 12.23387 

0.01608 10.49999 


XS 

19.77894 

17.95305 

15.00259 

12.28356 

10.59868 


YP XP 

-0.20109 19. 77817 

-0.16533 17.95241 

-0.09066 15.00232 

-0.00617 12.28394 

0.06570 10.59980 


XS 

19.72894 

17.94307 

15.00307 

12.33334 

10.69805 


YP XP 

-0.20070 19.72817 

-0.15804 17.94247 

-0.07250 15.00292 

0.02120 12.33392 

0.10389 10.69950 


XS 

19.67894 

17.93309 

15.00354 

12.38313 

10.79754 


YP XP 

-0.20031 19.67817 

-0.15226 17.93252 

-0.05866 15.00348 

0.04125 12.38387 

0.13126 10.79920 



z 


0* 90000 


ELEMENT 

YS 

XS 

YP 

XP 

1 

-0* 0^939 

19.62894 

-0.19988 

19.62817 

2 

-0.04781 

17.92310 

-0.14797 

17.92256 

3 

0.05119 

15.00401 

-0.04905 

15.00402 

4 

0.15402 

12.43299 

0.05418 

12.43382 

5 

0.24752 

10.89719 

0. 14814 

10.89899 


ELEMENT 

Z = 

1.00000 

YS 

XS 

YP 

XP 

1 


-0.09909 

19.57894 

-0. 19942 

19.57818 

2 


-0.04510 

17.91311 

-0.14510 

17.91258 

3 


0.05640 

15.00449 

-0.04363 

15.00453 

4 


0.16002 

12.48292 

0.06013 

12.48380 

5 


0.25468 

10.99705 

0.15469 

10.99891 



Z = 

1.10000 




ELEMENT 


YS 

XS 

YP 

XP 

1 


-0.09881 

19.52894 

-0.19894 

19.52818 

2 


-0.04386 

17.90311 

-0.14388 

17.90259 

3 


0.05762 

15.00499 

-0.04239 

15.00504 

4 


0.15948 

12.53293 

0.05914 

12.53380 

5 


0.25207 

11.09714 

0. 15098 

11.09897 


ELEMENT 

Z = 

1.20000 

YS 

XS 

YP 

XP 

1 


-0.09854 

19.47894 

-0.19844 

19.47818 

2 


-0.04411 

17.89311 

-0.14408 

17.89259 

3 


0.05484 

15.00550 

-0.04528 

15.00553 

4 


0.15232 

12.58302 

0.05122 

12.58384 

5 


0.23947 

11.19744 

0.13689 

11.19916 


ELEMENT 

Z = 

1.30000 

YS 

XS 

YP 

XP 

1 


-0.09826 

19.42894 

-0.19796 

19.42818 

2 


-0.04582 

17.88311 

-0.14577 

17.88257 

3 


0.04806 

15.00602 

-0.05235 

15.00601 

4 


0.13843 

12.63318 

0.03626 

12.63389 

5 


0.21659 

11.29792 

0.11208 

11.29944 


-a 

N3 


ELEMENT 

1 

2 

3 

4 

5 


FLEMENT 

1 

2 

3 

4 

5 


ELEMENT 

1 

2 

3 

4 

5 


ELEMENT 

1 

2 

3 

4 

5 


ELEMENT 

1 

2 

3 

4 

5 


Z = 1.40000 

YS XS 

-0.09796 19.37894 

-0.04901 17.87310 

0.03723 15.00655 

0.11764 12.68340 

0.18293 11.39853 


Z = 1.50000 

YS XS 

-0.09763 19.32894 

-0.05366 17.86309 

0.02234 15.00708 

0.08966 12.73363 

0.13783 11.49917 


Z = 1.60000 

YS XS 

-0.09728 19.27895 

-C. 05979 17.85307 

0.00327 15.00760 

0.05414 12.78383 

0.08034 11.59972 


Z = 1.70000 

YS XS 

-0.09690 19.22895 

-C. 06739 17.84304 

-0.02008 15.00808 

0.01059 12.83394 

0.00921 11.70000 


Z *= 1.80000 

YS XS 

-0.09650 19.17895 

-0.07647 17.83300 

-0.04779 15.00852 

-0.04163 12.88307 

-0.07736 11.79975 


YP 

-0.19750 

-0.14897 

-0.06363 

0.01404 

0.07601 


YP 

-0.19707 

-0.15367 

-0.07917 

-0.01574 

0.02794 


YP 

-0.19666 

-0.15988 

-0.09905 

-0.05351 

-0.03320 


-0.19627 
-0.16760 
-0.12342 
-0.09982 
-0. 10889 


YP 

-0.19592 

-0.17683 

-0.15239 

-0.15543 

-0.20119 


XP 

19.37819 

17.87255 

15.00646 

12.68393 

11.39975 


XP 

19.32819 

17.86251 

15.00689 

12.73393 

11.49997 


XP 

19.27819 

17.85245 

15.00727 

12.78383 

11.59995 


XP 

19.22819 

17.84238 

15.00759 

12.83355 

11.69949 


XP 

19.17819 

17.83229 

15.00782 

12.88300 

11.79828 



z 


1.90000 


ELEMENT 

YS 

xs 

1 

-0.09607 

19.12895 

2 

-0.08704 

17.82296 

3 

-0.08005 

15.00888 

4 

-0.10337 

12.93353 

5 

-0.18187 

11.89861 



Z = 

2.00000 


ELEMENT 


YS 

XS 

1 


-0.09561 

19.07895 

2 


-0.09910 

17.81289 

3 


-0.11702 

15.00914 

4 


-0.17575 

12.98275 

5 


-0.30808 

11.99604 



Z = 

2.10000 


ELEMENT 


YS 

XS 

1 


-0.09513 

19.02895 

2 


-0.11265 

17.80281 

3 


-0.15899 

15.00926 

4 


-0.26034 

13.03134 

5 


-0.46201 

12.09118 



Z = 

2.20000 


ELEMENT 


YS 

XS 

1 


-0.09463 

18.97896 

2 


-0. L2772 

L7. 79271 

3 


-0.20613 

15.00918 

4 


-0.35933 

13.07901 

5 


-0.65448 

12.18243 


ELEMENT 

YIS 

XI s 

ns 

YIP 

1 

-0.1025 

19.9997 

0.1584 

-0.2025 

2 

-0.1197 

18.0003 

0.1243 

-0.2186 

3 

-0.1591 

14.9991 

0.0600 

-0.2505 

4 

-0.2229 

11.9845 

0.0052 

-0.3033 

5 

-0.3460 

9.9690 

-0.0250 

-0.4185 




-a 

03 


YP XP 

-0.19558 19.12819 

-0.18757 17. 82218 

-0.18612 15.00794 

-0.22136 12.93205 

-0.31315 11.89588 


YP XP 

-0.19528 19.07819 

-0.19985 17.81205 

-0.22486 15.00791 

-0.29899 12.98050 

-0.44947 11.99158 


YP XP 

-0.19499 19.02819 

-0.21366 17.80189 

-0.26886 15.00769 

-0.39029 13.02810 

-0.61809 12.08420 


YP XP 

-0.19474 18.97819 

-0.22900 17.79169 

-0.31845 15.00722 

-0.49816 13.07445 

-0.83470 12.17141 


XIP ZIP YIC XIC ZIC 

19.9990 0.1584 -0.1525 19.9994 0.1584 

17.9979 0.1390 -0.1691 17.9992 0.1317 

14.9980 0.10Q7 -0.2048 14.9986 0.0804 

12.0096 0.0590 -0.2631 11.9971 0.0321 

10.0163 0.0250 -0.3822 9.9927 -0. 


ELEMENT 

YMS 

XMS 

ZMS 

YMP 

XMP 

ZMP 

YMC 

XMC 

ZMC 

I 

-0.09 8 7 

19.4997 

1.1584 

-0.1986 

19.4990 

1.1584 

-0.1487 

19.4994 

1.1584 

2 

-0.0438 

17.8999 

1.1317 

-0.1438 

17.8994 

1.1317 

-0.0938 

17.8998 

1.1317 

3 

0.0577 

15.0049 

1.0804 

-0.0423 

15.0049 

1.0804 

0.0077 

15.0050 

1.0804 

4 

0.1606 

12.4990 

1.0321 

0.0606 

12.4999 

1.0321 

0.1106 

12.4995 

1.0321 

5 

0.2547 

10.9971 

1.0000 

0.1547 

10.9989 

1.0000 

0.2047 

10. 9981 

1.0000 


ELEMENT 

YTS 

XTS 

ZTS 

YTP 

XTP 

ZTP 

YTC 

XTC 

ZTC 

1 

-0.0987 

19.4997 

1 • 1584 

-0.1986 

19.4990 

1.1584 

-0.1487 

19.4994 

1.1584 

2 

-0.0438 

17.8999 

1.1317 

-0.1438 

17.8994 

1.1317 

-0.0938 

17.8998 

1.1317 

3 

0.0577 

15.0049 

1.0804 

-0.0423 

15.0049 

1.0804 

0.0077 

15.0050 

1.0804 

4 

0.1606 

12.4990 

1.0321 

0.0606 

12.4999 

1.0321 

0.1106 

12.4995 

1.0321 

5 

0.2547 

10.9971 

1.0000 

0.1547 

10.9989 

1.0000 

0.2047 

10.9981 

1.0000 


REM ENT 

YOS 

XOS 

ZOS 

YOP 

XOP 

ZOP 

YOC 

XOC 

ZOC 

1 

-0.0948 

18.9998 

2.1584 

-0.1948 

18.9990 

2.1584 

-0.1448 

18.9994 

2.1584 

2 

-0.1184 

17.7989 

2.1390 

-0.2172 

17.7994 

2.1243 

-0.1678 

17.7992 

2.1317 

3 

-0.1593 

15.0093 

2.1007 

-0.2506 

15.0078 

2.0600 

-0.2049 

15.0086 

2.0804 

4 

-0.2241 

13.0115 

2.0590 

-0.3034 

12.98 30 

2.0052 

-0.2638 

12.9973 

2.0321 

5 

-0.34 3 7 

12.0201 

2.0250 

-0.4128 

11.9679 

1.9750 

-0.3782 

11.9940 

2.0000 


8LA0E SECTION COORDINATES (UNROTATEO) AT X = 11.0000 


ZIC 

ZMC 

ZTC 

ZOC 

ZIP 

0.0161 

l.OOOC 

1.0000 

1.9681 

0.0424 

YTC 

YMC 

YTC 

YOC 

YIP 

0.3126 

0.2046 

0.2046 

-0.4496 

-0.3519 

Tt 

TM 

TO 



0.09 47 

0.1000 

0.0743 


l 


- 0.0000 
0.1000 
0.2000 
0.3000 
0.4000 
0.5000 
0.6000 
0.7000 
0.8000 
0.9000 
1.0000 
1.1000 
1.2000 
1.3000 
1.4000 
1.5000 
1.6000 
1.7000 
1.8000 
1.9000 
2.0000 
2. lOOO 
2.2000 


ZMP 

ZTP 

ZOP 

ZIS 

ZMS 

1.0000 

1.0000 

1.9462 

-0.0102 

1.0001 

YMP 

YTP 

YOP 

YIS 

YMS 

0.1546 

0.1546 

-0.4796 

-0.2732 

0.2545 


YP 

YS 

-0.3977 

-0.2636 

-0.2938 

-0.1694 

-0.2040 

-0.0869 

-0.1263 

-0.0147 

-0.0591 

0.0483 

-0.0014 

0.1028 

0.0473 

0.1491 

0.0872 

0.1875 

0.1186 

0.2181 

0.1412 

0.2405 

0.1546 

0.2545 

0.1573 

0.2583 

0.1492 

0.2519 

0.1291 

0.2340 

0.0953 

0.2032 

0.045 ' 

0.1573 

-0.022V 

0.0936 

-0.1145 

0.0084 

-0.2349 

-0.1037 

-0.3925 

-0.2500 

-0.6003 

-0.4417 

-0.8798 

-0.6961 

-1.2713 

-1.0437 


ZTS ZOS ZCG 

1.0001 1.9900 1.0087 

YTS YDS YCG 

0.2545 -0.4196 0.0002 



2 IC 

0.04 79 
YIC 

-0.2407 

TI 

0.09 81 


ZMC 

1.0420 

YMC 

0.0889 

TM 

0.100C 


ZTC 

1.0420 

YTC 

0.0889 

TO 

0.0940 


BLADE SECTION COORDINATES (UNROTATED) AT X = 13.0000 


zoc 

ZIP 

ZMP 

ZTP 

ZOP 

ZIS 

ZMS 

2.0322 

0.0732 

1.0420 

1.0420 

2.0057 

0.0226 

1.0421 

YOC 

YIP 

YMP 

YTP 

YOP 

YIS 

YMS 

-0.2641 

-0.2827 

0.0390 

0.0390 

-0.3029 

-0.1987 

0.1390 


Z YP YS 


-0.0000 

-0.3358 

-0.2143 

0.1000 

-0.2645 

-3.1482 

0.2000 

-0.2024 

-0.0902 

0.3000 

-0.1484 

-0.0396 

0.4000 

-0.1020 

0.0041 

0.5000 

-0.0626 

0.0414 

0.6000 

-0.0299 

0.0725 

0.7000 

-0.0036 

0.0976 

0.8000 

0.0164 

0.1167 

0.9000 

0.0301 

0.1301 

1.0000 

0.0377 

0.1376 

1.1000 

0.0389 

0.1392 

1.2000 

0.0338 

0.1346 

1.3000 

0.0219 

0.1237 

1.4000 

0.0030 

0.1061 

1.5000 

-0.0234 

0.0815 

1.6000 

-0.0579 

0.0491 

1.7000 

-0.1013 

0.0084 

1.8000 

-0.1548 

-0.0416 

1.9000 

-0.2197 

-0.1022 

2.0000 

-0.2981 

-0.1751 

2.1000 

-0.3955 

-0.2655 

2.2000 

-0.5198 

-0.3796 


ZTS 

1.0421 

YTS 

0.1390 


ZOS 

2.0587 

YOS 

-0.2252 


ZCG 

1.0436 

YCG 

-0.0348 


-4 

U1 


ZIC 

0.0804 

YIC 

-0.2048 

TI 

0.1000 


ZMC 

1.0803 

YMC 

0.0079 

TM 

0.1000 


ZTC 

1.0803 

YTC 

0.0079 

TQ 

0.0999 


BLADE SECTION COORDINATES (UNROTATED) AT X = 15.0000 


ZOC 

ZIP 

ZMP 

ZTP 

ZOP 

ZIS 

ZMS 

2.0802 

0.1007 

1.0803 

1.0803 

2.0598 

0.0601 

1. 0803 

YOC 

YIP 

YMP 

YTP 

YOP 

YIS 

YMS 

0.2052 

-0.2504 

-0.0421 

-0.0421 

-0.2508 

-0.1591 

0.0579 


Z 

- 0.0000 

0.1000 

0.2000 

0.3000 

0.4000 

0.5000 

0.6000 

0.7000 

0.8000 

0.9000 

1.0000 

1.1000 

1.2000 

1.3000 

1.4000 

1.5000 

1.6000 

1.7000 

1.8000 

1.9000 

2.0000 

2.1000 

2.2000 


YP 

- 0.2980 

- 0.2508 

- 0.2089 

- 0.1722 

- 0.1403 

- 0.1132 

- 0.0906 

- 0.0724 

- 0.0585 

- 0.0489 

- 0.0435 

- 0.0422 

- 0.0451 

- 0 . 7)521 

- 0.0634 

- 0.0790 

- 0.0989 

-o;rm 

- 0.1524 

- 0.1862 

- 0.2251 

- 0.2692 

- 0.3189 


YS 

- 0.1867 

- 0.1417 

- 0.1018 

- 0.0667 

- 0.0363 

- 0.0103 

0.0114 

0.0288 

0.0421 

0.0513 

0.0566 

0.0578 

0.0550 

0.0483 

0.0374 

0.0225 

0.0034 

- 0.0200 

- 0.0477 

- 0.0801 

- 0.1172 

- 0.1593 

- 0.2067 


ZTS 

1.0803 

YTS 

0.0579 


ZOS 

2.1005 

YOS 

-0.1596 


ZCG 

1.0785 

YCG 

-0.0697 



-a 

o> 


7 TC ZMC 

0.1317 1.1334 

YTC YMC 

-0.1691 -0.0972 

T I TM 

0.0999 0.1000 


ZTC 

1.1334 

YTC 

-0.0972 

TO 

0.1000 


BLADE SECTION COORDINATES (UNROTATED) AT X = 18.0000 


zoc 

ZIP 

ZMP 

ZTP 

ZOP 

ZIS 

ZMS 

ZTS 

2.1360 

0.1391 

1.1334 

1.1334 

2.1298 

0.1243 

1.1334 

1.1334 

YOC 

YIP 

YMP 

YTP 

YOP 

YIS 

YMS 

YTS 

-0.1644 

-0.2185 

-0.1472 

-0.1472 

-0.2140 

-0.1197 

-0.0472 

-0.0472 


Z YP YS 


-0.0000 

-0.2412 

-0.1399 

0.1000 

-0.2246 

-0.1234 

0.2000 

-0.2096 

-0.1087 

0.3000 

-0.1964 

-0.0955 

0.4000 

-0.1848 

-0.0841 

0.5000 

-0.1748 

-0.0742 

0.6000 

-0.1664 

-0.0659 

0.7000 

-0.1595 

-0.0592 

0.8000 

-0.1542 

-0.0540 

0.9000 

-0.1504 

-0.0502 

1.0000 

-0.1481 

-0.0480 

1.1000 

-0.1472 

-0.0472 

1.2000 

-0.1478 

-0.0478 

1.3000 

-0.1497 

-0.0497 

1.4000 

-0.1530 

-0.0530 

1.5000 

-0.1576 

-0.0576 

1.6000 

-0.1635 

-0.0635 

1.7000 

-0.1706 

-0.0705 

1.8000 

-0.1789 

-0.0787 

1.9000 

-0.1884 

-0.0880 

2.0000 

-0.1989 

-0.0984 

2.1000 

-0.2104 

-0.1097 

2.2000 

-0.2228 

-0.1219 


ZDS 

2.1422 

YOS 

-0.1148 


2 IC 

0.15 23 
YIC 

-0.1565 

TI 

0.0999 


ZMC ZTC 

1.1584 1.1584 

YMC YTC 

-0.1487 -0.1487 

TM TO 

0.1000 0.0992 


BLADE SECTION COORDINATES (UNROTATED! AT X = 19.5000 


znc 

ZIP 

ZMP 

ZTP 

ZOP 

ZIS 

ZMS 

2.1695 

0.1540 

1.1584 

1.1584 

2.1726 

0.1505 

1.1584 

YOC 

YIP 

YMP 

YTP 

YOP 

YIS 

YMS 

-0.1326 

-0.2064 

-0.1987 

-0.1987 

-0.1821 

-0.1066 

-0.0987 


Z YP YS 


0.0000 

-0.2134 

-0.1138 

0.1000 

-0.2087 

-0.1088 

0.2000 

-0.2047 

-0.1046 

0.3000 

-0.2015 

-0.1012 

0.4000 

-0.1990 

-0.0986 

0.5000 

-0.1973 

-0.0968 

0.6000 

-0.1962 

-0.0956 

0.7000 

-0.1956 

-0.0951 

0.8000 

-0.1957 

-0.0951 

0.9000 

-0.1961 

-0.0957 

1.0000 

-0.1969 

-0.0966 

1.1000 

-0.1980 

-0.0978 

1.2000 

-0.1992 

-0.0993 

1.3000 

-0.2004 

-0.1007 

1.4000 

-0.2014 

-0.1019 

1.5000 

-0.2022 

-0.1028 

1.6000 

-0.2024 

-0.1031 

1.7000 

-0.2019 

-0.1027 

1.8000 

-0.2005 

-0.1013 

1.9000 

-0.1979 

-0.0986 

2.0000 

-0.1938 

-0.0944 

2.1000 

-0.1878 

-0.0884 

2.2000 

-0.1797 

-0.0801 


ZTS 

1.1584 

YTS 

-0.0987 


ZOS 

2.1664 

YOS 

-0.0831 


ZCG 

1.1309 

YCG 

- 0.1220 


ZCG 

1.1571 

YCG 

-0.1482 



BLADE SECTION COORDINATES (ROTATED) AT X = 11.0000 


GAMMA 

T I 

MSP) 

L-BAR 

AREA 

IMTN ILLCG 


PHLCG 

I(LL) 

PHL 



3.7191 

0.0947 

1.0177 

1.0143 

0.2252 

0.8117 E-02 0.8394E 

o 

(V 

0 

1 

. 4536E-02 

0.4830E-01 

0.1007 



TM 

TO 

H ( SP) 

H-BAR 

BETA 

IMAX IHHCG 



I(HH) 




0.1000 

0.0743 

0.4236 

0.4209 

3.4890 

0 . 8279E-0 1 0.8251F 

-01 


0.3142 




L C 1 C ) 

KMC) 

L ( TC ) 

L ( OC ) 

KIP) 

UMP) L(TP) 

L ( OP ) 

L ( IS) 

L ( MS ) 

KTS) 

L ( OS) 

L(CG) 

0.0473 

0.9956 

0.9956 

2.0041 

0.0761 

0.9989 0.9989 1 

.9842 

0.0185 

0.9924 

0.9924 

2.0240 

1.0176 

H(IC) 

HI MC ) 

H ( TC ) 

H ( OC ) 

H ( I P) 

H(MP) H ( T P ) 

H( OP ) 

H( IS) 

H{ MS) 

H ( TS ) 

H ( OS) 

H ( CG) 

0.0473 

0.6272 

0.6272 

0.0371 

0.0098 

0.5773 0.5773 0 

.0058 

0.0849 

0.6770 

0.6770 

0.0685 

0.4238 


L HP HS 


0 . 

0.0473 

0.0473 

0.0473 

-0.0254 

0.1168 

0.0500 

-0.0221 

0.1198 

0.1000 

0.0378 

0.1734 

0. 1500 

0.0932 

0.2233 

0.2000 

0.1444 

0.2699 

0.2500 

0.1919 

0.3134 

0.3000 

0.2360 

0.3541 

0.3500 

0.2770 

0.3921 

0.4000 

0.3150 

0.4275 

0.4500 

0.3503 

0.4605 

0.5000 

0.3830 

0.4912 

0.5500 

0*4131 

0.5196 

0.6000 

0.4408 

0.5459 

0.6500 

0.4660 

0.5699 

0.7000 

0.4890 

0.5919 

0.7500 

0.5096 

0.6117 

0.8000 

0.5279 

0.6293 

0.8500 

0.5439 

0.6449 

0.9000 

0.5575 

0.6582 

0.9500 

0.5689 

0.6695 

1.0000 

0.5775 

0.6781 

1.0500 

0.5834 

0.6842 

1.1000 

0.5866 

0.6877 

1.1500 

0.5872 

0.6887 

1.2000 

0.5848 

0.6869 

1.2500 

0.5797 

0.6824 

1,3000 

0.5713 

0.6748 

1.3500 

0.5599 

0.6642 

1.4000 

0.5447 

0.6502 

1.4500 

0.5262 

0.6328 

1.5000 

0.5035 

0.6114 

1.5500 

0.4769 

0.5862 

1.6000 

0.4454' 

0.5564 

1.6500 

0.4094 

0.5223 

1.7000 

0.3677 

0.4827 

1.7500 

0.3207 

0.4381 

1.8000 

0.2670 

0.3869 

1.8500 

0.2068 

0.3298 

1.9000 

0.1387 

0.2649 

1.9500 

0.0625 

0.1926 

2.0000 

-0.0230 

0.1112 

2.0041 

-0.0305 

0.1041 

2.0412 

0.0371 

0.0371 



BLADE SECTION COORDINATES (ROTATED) AT X = 13.0000 


GAM M^ 

T I 

L ( SP) 

L-BAR 

AREA 

imtn illcg 

PHLCG 

I (LL) 

PHL 



0.61 50 

0.0981 

i .0074 

1.0433 

0.2187 

0 • 2804F-02 0.2816P-0? 0 

• 957 IF- 03 

0. 1845E-01 

0.6196E- 

-01 


TM 

TO 

H( SP) 

H-BAR 

BETA 

IMAX IHHCG 


HHH) 




n.iooo 

0.0940 

0.3001 

0.2673 

0.706? 

0.R047E-01 0.8046F-01 


0.3185 




L ( I C ) 

L ( M C ) 

l ( TC ) 

L ( OC ) 

L( IP) 

L ( MP ) L ( T P ) L(HP) 

L ( IS) 

L(MS) 

L ( TS ) 

L(OS) 

L ( CG) 

0 .0491 

1 .0396 

1.0396 

2.0335 

0.0748 

1.0401 1.0401 2.0074 

0.0233 

1.0391 

1.0391 

'2.0596 

1.0425 

H( If. ) 

H ( MC ) 

H ( TC ) 

H { OC ) 

H ( I ° ) 

H( MP) H(TP> H(DP) 

H( IS) 

H(MS) 

H ( T S ) 

H ( QS ) 

H ( CG) 

0.0491 

0.3893 

0.3893 

0.0470 

0.0073 

0.3394 0.3394 0.0079 

0.0908 

0.4394 

0.4394 

0.0861 

0.2657 


L HP HS 


0 . 

0.0491 

0.0491 

0.0491 

-0.0112 

0.1086 

0.0500 

-0.0106 

0.1092 

0.1000 

0.0248 

0.1420 

0. 1500 

0.0577 

0.1726 

0.2000 

0.0885 

0.2014 

0.2500 

0.) 172 

0.2282 

0.3000 

0.1439 

0.2533 

0.3500 

0.1687 

0.2766 

0.4000 

0.1917 

0.2983 

0.4500 

0.2128 

0.3183 

0.5000 

0.2323 

0.3367 

0.5500 

0.2500 

0.3535 

0.6000 

0.2661 

0.3688 

0.6500 

0.2806 

0.3826 

0.7000 

0.2935 

0.3949 

0.7500 

0.3048 

0.4057 

0.8000 

0.3145 

0.4150 

0.8500 

0.3226 

0.4229 

0.9000 

0.3292 

0.4293 

0.9500 

0.3343 

0.4342 

1.0000 

0.3378 

0.4377 

1.0500 

0.3397 

0.4397 

1.1000 

0.3400 

0.4402 

1.1500 

0.3387 

0.4391 

1.2000 

0.3357 

0.4365 

1.2500 

0.3312 

0.4323 

1.3000 

0.3248 

0.4265 

1.3500 

0.3168 

0.4191 

1.4000 

0.3070 

0.4099 

1.4500 

0.2953 

0.3989 

1.5000 

0.2816 

0.3861 

1.5500 

0.2660 

0.3715 

1.6000 

0.2482 

0.3548 

1.6500 

0.2283 

0.3362 

1.7000 

0.2060 

0.3152 

1.7500 

0.1813 

0.2921 

1.8000 

0.1539 

0.2664 

1.8500 

0.1238 

0.2382 

1.9000 

0.0905 

0.2071 

1.9500 

0.0546 

0.1736 

2.0000 

0.0141 

0.1359 

2.0335 

-0.0148 

0.1089 

2.0805 

0.0470 

0.0470 



SLADE SECTION COORDINATES (ROTATED) AT X = 15.0000 


GAMMA 

T I 

L( SP) 

L-BAR 

AREA 

1 MTN TLLCG 

PHLCG 

I (LL ) 

PHL 



0.0103 

0.1000 

0.9784 

1.0487 

0.2133 

0.1 176E-02 0 . 1 176E-02 0 

.5622E-04 

0 . 85 8 5E-02 

0.4175E- 

-01 


TM 

Tn 

H ( SP ) 

H-BAR 

BETA 

I MAX IHHCG 


I (HH ) 




0.1000 

0.0999 

0.2549 

0.1864 

0.4202F- 

-01 0.7782E-01 O.7702E-O1 


0.3124 




L ( T C ) 

L(MC) 

L ( TC ) 

L(OC) 

L ( I P ) 

L ( MP ) L ( TP ) HOP) 

L ( IS) 

L ( MS) 

L(TS) 

L ( OS ) 

L ( CG) 

0.0500 

1.0499 

1.0499 

2.0498 

0.0703 

1.0499 1.0499 2.0294 

0.0296 

1.0498 

1 .0498 

2.0701 

1.0481 

H(1C) 

H ( MC ) 

H ( TC ) 

H( no 

HUP) 

H ( MP) H(TP) H(OP) 

H( IS) 

H ( MS ) 

H ( TS ) 

H ( OS ) 

H ( CG ) 

0.05 00 

0.2628 

0.2628 

0.0499 

0.0043 

0.2128 0.2128 0.0043 

0.0957 

0.3128 

0.3128 

0.0956 

0.1853 


L HP HS 


0 . 

0.0500 

0.0500 

0.0500 

-0.0048 

0.1046 

0.0500 

-0.0048 

0.1046 

0.1000 

0.0173 

0.1257 

0.1500 

0.0381 

0.1455 

0.2000 

0.0576 

0.1641 

0.2500 

0.0758 

0.1816 

0.3000 

0.0928 

0.1978 

0.3500 

0.1086 

0.2129 

0.4000 

0.1232 

0.2269 

0.4500 

0.1367 

0.2398 

0.5000 

0.1490 

0.2516 

0.5500 

0.16*02 

0.2623 

0.6000 

0.1703 

0.2720 

0.6500 

0.1792 

0.2806 

0.7000 

0.1871 

0.2882 

0.7500 

0.1940 

0.2947 

0.8000 

0.1997 

0.3002 

0.8500 

0.2044 

0.3048 

0.9000 

0.2081 

0.3083 

0.9500 

0.2107 

0.3108 

1.0000 

0.2123 

0.3123 

1.0500 

0.2128 

0.3128 

1.1000 

0.2123 

0.3123 

1.1500 

0.2108 

0.3109 

1.2000 

0.2082 

0.3084 

1.2500 

0.2045 

0.3049 

1.3000 

0.1999 

0.3004 

1.3500 

0.1941 

0.2949 

1.4000 

0.1873 

0.2884 

1.4500 

0.1794 

0.2808 

1.5000 

0.1705 

0.2722 

1.5500 

0.1604 

0.2625 

1.6000 

0.1492 

0.2518 

1.6500 

0.1369 

0.2400 

1.7000 

0.1234 

0.2271 

1.7500 

0.1088 

0.2131 

1.8000 

0.0930 

0.1980 

1.8500 

0.0759 

0.1817 

1.9000 

0.0576 

0.1642 

1.9500 

0.0381 

0.1456 

2.0000 

0.0172 

0.1256 

2.0498 

-0.0049 

0.1046 

2.0997 

0.0499 

0.0499 


CD 



BLADE SECTION COORDINATES (ROTATED) AT X = 18.0000 


GAMMA 

T I 

KSP) 

L-BAR 

AREA 

I M IN ILLCG 

PHLCG 

I (LL ) 

PHL 



0.1340 

0.0999 

0.9275 

1. 0494 

0.2076 

0. 2755 E-03 0.2758E-03 -0 

• 1358E-03 

0 .21 36E-Q2 

0.2049E- 

•01 


TM 

TO 

H ( SP ) 

H-BAR 

BETA 

IMAX IHHCG 


I ( HH ) 




0.10 00 

0. 1000 

0.2170 

0.0946 

-0.1052 

0.7426E-01 0 . 7 42 6 E- 0 1 


0.3029 




LUC) 

L ( MC ) 

L(TC) 

L ( OC ) 

LIIPI 

L(MP) LCTP) L ( OP ) 

L ( IS) 

L ( MS) 

L ( TS) 

L ( OS) 

L(CG) 

0.05 00 

1.0518 

1 *051 B 

2.0542 

0.0572 

1.0517 1.0517 2.0479 

0.0427 

1.0519 

1.0519 

2.0606 

1.0493 

H( IC) 

H ( MC ) 

H ( TC ) 

H( OC) 

HUP) 

H ( MP ) HUP) H( np ) 

H( IS) 

H ( MS) 

HUS) 

H ( OS) 

H(CG) 

0.05 00 

0.1195 

0.1195 

0.0500 

0.0005 

0.0695 0.0695 0.0004 

0.0994 

0.1695 

0.1695 

0.0996 

0.0947 


L HP HS 


0. 

0.0500 

0.0500 

0.0500 

-0.0006 

0. 1005 

0.0500 

-0.0006 

0.1005 

0. 1000 

0.0067 

0.1077 

0.1500 

0.0136 

0.1145 

0.2000 

0.0201 

0.1209 

0.2500 

0.0261 

0.1269 

0.3000 

0.0317 

0.1324 

0.3500 

0.0370 

0.1376 

0.4000 

0.0418 

0.1424 

0.4500 

0.0462 

0.1467 

0.5000 

0.0502 

0.1507 

0.5500 

0.0539 

0.1543 

0.6000 

0.0571 

0.1575 

0.6500 

0.0600 

0.1603 

0.7000 

0.0625 

0.1628 

0.7500 

0.0646 

0.1648 

0.8000 

0.0663 

0.1665 

0.8500 

0.0677 

0.1678 

0.9000 

0.0687 

0.1688 

0.9500 

0.0693 

0.1694 

1.0000 

0.0696 

0.1696 

1.0500 

0.0695 
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